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Technical  findings: 

1.  Bonding  effects  on  the  slip  differences  in  the  B1  monocarbides 
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Thompson1* 
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Street,  Philadelphia,  PA  19104 
*  Correspondence  to:  gthompson@eng.ua.edu 
** Accepted  in  PRL,  2015 

Abstract: 

Differences  in  plasticity  are  usually  attributed  to  significant  changes  in  crystalline  symmetry  or 
the  strength  the  inter-atomic  bonds.  In  the  B1  monocarbides,  differences  in  slip  planes  exist  at 
low  temperatures  despite  having  the  same  structure  and  very  similar  bonding  characteristics. 

Our  experimental  results  demonstrate  concretely  that  HfC  slips  on  {110}  planes  while  TaC  slips 
on  {111}  planes.  Density  functional  theory  calculations  rationalize  this  differences  through  the 
formation  of  an  intrinsic  stacking  fault  on  the  {111}  planes,  formation  of  Shockley  partials,  and 
enhanced  metallic  bonding  because  of  the  valence  filling  of  electrons  between  these  transitional 
metal  carbides. 

Main  body: 

Plasticity  and  ductility  are  generally  associated  with  crystalline  symmetry  and  the 
strength  of  the  inter-atomic  bonds  [1].  For  example,  face  centered  cubic  (FCC)  and  hexagonal 
close  packed  (HCP)  materials  have  equivalent  packing  with  similar  types  of  bonding  but  show 
dramatic  differences  in  deformation  because  of  the  availability  of  slip  systems  associated  with 
their  symmetry.  Alternatively,  dramatic  differences  in  bonding,  such  as  those  between  ceramics 
and  metals,  can  also  have  a  similar  impact  on  dislocation  slip  and  ductility,  but  such  differences 
are  confounded  by  variances  in  symmetry.  What  can  be  even  more  perplexing  are  materials  that 
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exhibit  equivalent  symmetry  and  bonding  but  show  differences  in  defonnation  responses.  In  this 
letter,  we  explore  how  a  subtle  difference  in  bonding  behavior  regulates  the  slip  response  in  two 
equivalently  structured  and  similarly  bonded  monocarbides  by  combining  density  functional 
theory  (DFT)  and  experimental  results.  Our  findings  debunk  how  prior  hard  sphere  based  models 
are  not  able  to  rationalize  experimental  findings.  Through  the  use  of  these  modem  computational 
tools,  we  elucidate  how  bonding,  along  with  symmetry,  regulates  slip  in  this  class  of  transition 
metal  monocarbides  (TMMCs). 

The  group  IVB  and  VB  monocarbides  provide  a  perfect  test-bed  to  explore  subtle 
differences  in  bonding  and  its  role  in  mechanical  properties.  Both  groups  of  TMMCs  fonn  the 
Bl,  or  rocksalt  structure,  which  is  an  array  of  FCC  metal  atoms  with  the  carbon  atoms  filling 
every  octahedral  interstice  [2,3].  The  bonding  in  these  materials  is  a  mix  of  covalent,  metallic 
and  ionic  resulting  in  high  hardness  as  well  as  good  thermal  and  electrical  conductivity  [4-7]. 

The  materials  are  generally  brittle  at  room  temperature  and  exhibit  remarkable  ductility  at 
elevated  temperatures  [8-10]. 

The  last  comprehensive  summary  of  plastic  deformation  of  TMMCs  was  written  by 
Rowcliffe  in  1984  [11],  where  most  slip  systems  were  detennined  using  hardness  anisotropy. 

This  work  reported  that  the  dominant  slip  systems  in  the  group  IVB  monocarbides  are 
a/2<l  1 0>  {110}  while  the  group  VB  monocarbides  revealed  a/2<l  1 0>  {111}  at  room 
temperature.  Rowcliffe  et  al.  [12]  also  noted  that  the  group  IVB  TMMCs  were  significantly  more 
brittle  than  the  group  VB  TMMCs,  which  is  attributed  to  the  number  of  available  independent 
slip  systems  in  the  dominant  family,  i.e.  two  for  the  a/2<l  10>{  1 10}  and  five  for  a/2<l  10>{  1 11}. 
The  supplemental  material  of  this  letter  provides  a  review  and  summary  table  which  shows  both 
slip  systems  are  found  in  all  of  the  materials  [13],  However,  it  has  been  established  from  past 
experiments  [8,1 1,14-16]  that  one  family  of  slip  system  appears  to  be  dominant  for  each  type  of 
carbide,  i.e.  the  {110}  in  the  group  IVB  TMMCs  and  {111}  in  the  group  VB  TMMCs.  Though 
these  differences  have  been  known  for  decades,  the  change  in  slip  systems  across  the  TMMCs 
has  perplexed  scientists  for  equally  as  long. 

To  eliminate  differences  in  homologous  temperatures,  we  pick  the  two  monocarbides  with 
nearly  equivalent  melting  temperatures,  HfC  (-3900  °C)  and  TaC  (-3880  °C),  to  conduct  our 
studies.  The  deformation  of  TaC  has  been  well  investigated  across  a  number  of  temperatures 
establishing  its  preference  to  slip  on  {111}  planes  through  indentation  experiments  and  elevated 
temperature  creep  studies  [11,16].  Despite  the  well-known  preference  for  slip  on  {111}  planes, 
the  mechanism  controlling  slip  is  not  well  known.  Authors  have  suggested  a  simple  high  lattice 
friction  model  with  dislocations  splitting  into  Shockley  partials,  a  synchro-shear  mechanism 
regulated  by  carbon  diffusion,  as  well  as  motion  via  zonal  dislocations  to  rationalize 
experimental  observations  [9,1 1,17];  the  latter  two  being  unlikely  to  occur  during  room 
temperature  defonnation.  In  addition,  as  the  carbon  concentration  in  TaCi-x  is  reduced,  hardness 
anisotropy  measurements  indicated  that  {110}  slip  becomes  more  prevalent  [18],  In  contrast, 
there  has  been  limited  studies  in  HfC  and  few  theories  that  explain  the  {110}  slip  in  this  system 
or  other  similar  group  IVB  carbides.  The  typical  explanation  has  been  that  these  materials  have 
more  covalent  bonding  preventing  {111}  slip  [11,19].  As  the  temperature  increases,  one  can 
expect  that  multiple  slip  systems  would  become  more  favorable  and  a  dominate  slip  system,  as 
noted  in  lower  homologous  temperatures,  may  not  be  as  apparent. 
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To  verify  these  results  obtained  decades  ago,  we  perfonned  microindentation  tests,  Fig. 
1(a),  in  nearly  stoichiometric  HfC  and  TaC  at  room  temperature  with  the  corresponding 
dislocations  observed  in  the  TEM  micrographs  of  Figs.  1(b)  and  1(c).  The  TEM  results  confirm 
large  amounts  of  dislocation  plasticity  within  the  TaC  grains  under  the  indents  whereas  the 
plasticity  in  HfC  was  much  more  limited  and  confined  to  large  grains  directly  under  the  indents 
with  estimated  dislocation  densities  being  -1.10  x  1014  m"2  (TaC)  and  -8.83  x  10Ll  nT2  (HfC) 
under  the  indents  and  -1.32  xl014(TaC)  and  -3.12  xlO11  m"2  (HfC)  away  from  the  indents;  the 
experimental  details  can  be  found  in  the  supplementary  infonnation  (SI)  [13]. 

Dynamical  diffraction  TEM  analysis  confirmed  the  slip  system  as  a/2<l  10>{  1 1 1}  in  TaC 
and  a/2<l  10>{  1 10}  in  HfC  [13].  LECO  analysis  confirmed  nearly  stoichiometric  and  equal 
amounts  of  carbon  in  the  samples,  eliminating  this  as  a  cause  for  the  change  in  slip  planes. 

These  results,  in  addition  to  the  previous  observations  in  single  crystals  [12],  suggest  the  choice 
of  slip  planes  and  the  limited  ductility  in  HfC  compared  to  TaC  are  intrinsic  properties. 


HfC 

200  nm 

50  nm 


FIG.  1.  (a)  Representative  schematic  and  optical  micrograph  of  the  10x10  array  of  Vickers  indents 
in  a  carbide  specimen.  TEM  bright  field  micrographs  of  typical  dislocation  structures  found  in  HfC 
(b)  and  TaC  (c). 


It  is  generally  thought  that  slip  will  occur  in  the  closes  packed  directions  on  the  most 
widely  spaced  planes  This  suggests  slip  on  {111}  planes  in  FCC  metals  and  {110}  slip  in  BCC 
metals  and  is  thought  to  control  the  competition  between  basal  and  prims  slip  in  HCP  metals. 
The  other  factor  complicating  slip  is  the  existence  of  fault  planes  [20].  However,  this  can  be 
complicated  in  materials  with  mixed  bonding.  The  predominately  held  belief  with  regard  to  the 
differences  in  slip  for  group  IVB  and  VB  carbides  was  based  on  the  hard  sphere  model.  First 
proposed  to  be  viable  in  ionic  and  covalent  materials  by  Van  Der  Walt  and  Sole  [21],  Rowcliffe 
and  Hollox  [18]  supported  its  applicability  to  carbides.  The  hard  sphere  model  relates  preferred 
slip  to  the  interatomic  spacing  within  the  structure  as  determined  by  the  ratio  of  the  radii  of  the 
carbon  and  metal  atoms.  It  states  that  perfect  slip  on  the  {110}  plane  is  preferred  in  covalently 
bonded  materials  with  a  radius  ratio  (r/R)  <  0.414,  on  the  {111}  plane  when  0.414  <  r/R  < 
0.633,  and  on  the  {100}  when  r/R  >  0.633.  It  also  only  predicts  partial  slip  on  the  {111}  plane  at 
r/R  >  0.732.  The  radius  ratios  of  all  the  group  IVB  and  VB  carbides,  as  reported  by  Toth  [4], 
along  with  the  hard  sphere  model  slip  predictions  and  the  experimentally  reported  and  observed 
active  slip  systems  are  tabulated  in  the  SI  for  the  reader.  These  results  suggest  that  all  of  the 
TMMCs  should  deform  via  {111}  slip  despite  the  observations  of  the  {110}  slip  dominance  in 
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the  group  IVB  TMMCs.  Given  the  close  radius  ratio  between  TiC  and  the  group  VB  TMMCs, 
the  model  has  difficulty  even  predicting  trends.  It  is  evident  that  the  hard  sphere  model  does  not 
adequately  explain  or  predict  favorable  slip  within  the  TMMCs  even  though  it  has  been 
propagated  for  several  years  in  contrast  to  prior  experimental  data.  Moreover  it  further  breaks 
down  when  partial  slip  is  added  to  the  discussion  as  will  be  elucidated  shortly. 

To  provide  for  a  more  accurate  understanding  of  slip,  we  have  employed  first  principle 
DFT  calculations  of  the  generalized  stacking  fault  (GSF)  energies  (see  SI).  The  GSF  energy 
surfaces  for  perfect  slip  for  {110}  and  {111}  planes  are  plotted  in  Fig.  2(a).  It  was  found  that  in 
either  TMMC  system,  {110}  slip  was  more  energetically  favorable  than  {111}  when  compared 
along  the  <11 0>  direction.  Moreover  TaC  was  shown  to  have  a  lower  GSF  than  HfC,  which 
helps  explain  the  higher  dislocation  densities  noted  in  the  material  between  the  two  carbides,  Fig. 
1 .  Though  TaC  has  a  lower  GSF  energy  than  HfC,  this  does  not  explain  its  dominance  on  {111} 
planes,  as  perfect  slip  on  these  planes  are  nearly  two-fold  higher  in  energy  than  {110}  perfect 
slip,  Fig.  2(a).  The  choice  of  slip  planes  rather  lies  with  the  presence  of  an  intrinsic  fault,  Fig. 
2(b),  and  ability  of  the  perfect  dislocations  to  split  into  Shockley  partials,  Figs.  2(c)  and  2(d). 
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Fractional  shear  along  [110]  Fractional  shear  along  [1 12](1 11) 


FIG.  2.  DFT  calculated  generalized  stacking  fault  energy  curves  for  fractional  shear  along  (a)  [110] 
on  both  (111)  and  (110)  planes  and  (b)  for  the  partial  dissociation  [1 12](1 1 1)  for  both  HfC  and 
TaC.  Note  the  intrinsic  stacking  fault  in  TaC  which  is  absent  in  HfC.  Associated  PN  model  misfit 
density  plots  for  (c)  HfC  and  (d)  TaC,  with  the  distance  normalized  by  the  Burgers  vector.  The 
preference  for  splitting  into  partials  is  evident  in  TaC  and  not  HfC. 
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It  is  well  known  that  dislocations  in  close-packed  metals,  amongst  others,  dissociate  into 
partial  dislocations  to  reduce  the  elastic  energy  of  the  system  [20].  To  explore  this  possibility  in 
the  TMMCs,  we  computed  the  <1 12>{  1 11}  GSF  curves  for  HfC  and  TaC,  Fig.  2(b).  In  this 
configuration  a  local  ID  minima  in  TaC,  known  as  an  intrinsic  stacking  fault  (ISF),  is  noted  but 
is  absent  in  HfC.  This  ISF  was  noted  to  be  present  in  all  of  the  group  VB  TMMCs  but  absent  in 
the  group  IVB  TMMCs  (see  SI)  [13],  indicting  unifonnity  across  the  two  family  classes  of 
monocarbides  [13].  The  ISF,  when  present,  represents  a  metastable  minimum  energy 
configuration  in  the  energy  surface,  which  may  result  in  the  splitting  of  perfect  dislocations  into 
partial  dislocations.  The  unstable  stacking  fault  energy  (USF)  is  the  maximum  energy  along  the 
minimum  energy  path  and  can  indicate  the  relative  ductility  of  a  material  [22,23].  In  the  hard 
sphere  model,  an  ISF  is  only  present  for  radius  ratios  or  r/R  >  0.732,  when  the  hard  sphere  model 
predicts  {001}  slip.  The  presence  of  an  ISF  in  TaC,  which  would  correspond  to  {001}  slip  in  the 
hard  sphere  model,  is  a  direct  indication  of  a  fundamental  breakdown  of  the  hard  sphere  model. 
Moreover,  the  presence  of  a  stable  ISF  also  likely  contributes  to  the  phase  stability  of  faulted- 
phases,  such  as  Ta4C3,  in  sub-stoichiometric  tantalum  carbides  but  such  phases  are  not  observed 
in  the  hafnium  carbides  [24]. 

These  DFT  calculations  shed  new  light  into  the  differences  in  the  room  temperature  slip 
behavior  between  the  two  carbides  that  has  been  lacking  for  years.  Perfect  slip  on  the  {110} 
surface  is  energetically  most  favorable  for  HfC,  which  is  in  agreement  with  the  results  reported 
here  and  elsewhere  [11].  In  contrast,  TaC  can  supersede  the  {110}  slip  because  of  its  favorable 
ISF  on  the  {111}  plane,  which  is  absent  in  HfC. 

To  confirm  that  dislocations  dissociate  into  partials  in  TaC  on  the  {111}  plane,  we  used  a 
Peierls-Nabarro  (PN)  model  to  compute  the  potential  for  dislocation  dissociation.  These  results 
are  shown  in  Figs.  2(c)  and  2(d)  and  confirm  our  expectations  based  on  the  GSF  curves  alone.  A 
TaC  dislocation  will  split  slightly  into  partial  dislocations  on  the  {111}  plane  as  shown  by  the 
separate  peaks  in  Fig.  2(d).  In  contrast,  the  lack  of  the  ISF  prevents  dislocations  from  splitting  in 
HfC,  also  confirmed  by  our  PN  model  in  Fig.  2(d).  Hence,  the  dissociation  of  the  dislocations 
helps  stabilize  them  on  the  {111}  plane  and  generally  lowers  their  barrier  to  motion  making 
{111}  slip  more  prevalent  in  TaC  than  HfC. 

To  fully  understand  how  the  DFT  calculations  can  provide  an  alternate  explanation  to  the 
hard  sphere  model  concepts,  one  must  consider  the  nature  and  directionality  of  the  bonding 
within  the  carbides.  Fig.  3  shows  the  isoconcentration  charge  surfaces  computed  from  our  DFT 
simulations,  providing  a  visual  aid  in  highlighting  the  differences  in  the  bonding  between  the 
carbides  during  the  shear  process.  In  all  renderings,  the  isocharge  surfaces  were  equivalent  for 
direct  visual  comparison  purposes  [13];  see  SI  for  details 
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FIG.  3.  DFT  isoconcentration  charge  surfaces  for  HfC  and  TaC  in  a  zero  and  fraction  shear  state 
for  various  slip  conditions.  The  dash  arrow  indicates  the  fault  plane  highlighted  by  the  solid  box 
region.  The  top  row  are  the  unsheared  conditions  with  the  bottom  row  are  the  designated  shear 
indicated  by  the  fraction  value.  The  compass  arrows  indicate  the  viewing  prospective. 


In  the  un-sheared  state,  we  can  clearly  see  a  difference  the  bonding  that  occurs  in  these 
two  carbides.  Notably,  the  isocharge  surfaces  are  more  diffuse  between  the  atoms  in  TaC  than  in 
HfC  regardless  of  the  isocharge  value  used  in  the  visualization.  This  indicates  a  less  directional 
bond  and  thus  more  metallic-like  character.  This  is  a  direct  consequence  of  the  extra  valence 
electron  contributed  by  tantalum,  a  group  VB  transition  metal  atom,  to  the  Ta-C  bond.  The 
difference  in  bonding  is  further  highlighted  as  the  crystal  shears,  with  the  bonding  behaving  in  a 
more  localized  fashion  in  HfC  than  TaC,  i.e.  compare  the  charge  distribution  between  the  two 
shear  planes  (dashed  arrows)  for  the  0 . 5  shear  along  the  <1 10>{  110}. 

The  less  directional  nature  of  the  bonds  in  TaC  also  allows  for  the  stabilization  of  the 
ISF.  The  isocharge  surfaces  at  a  fractional  shear  of  0.35  in  the  a/6<l  12>  direction  on  the  {111} 
planes,  Fig.  3,  exhibit  a  rotation  of  the  bonds  for  both  TaC  and  HfC.  As  the  top  half  of  the 
crystal  shears  relative  to  the  bottom  half  via  a  Shockley  partial  dislocation,  the  bonds  between 
the  metal  atoms  and  carbon  atoms  in  this  plane  undergo  a  60  degree  rotation  about  the  plane 
normal.  This  bonding  rotation  increases  the  bond  energy  which  will  be  a  function  of  the  angular 
nature  of  the  bonds  themselves.  The  less  directional  nature  of  the  Ta-C  bonds  helps  to  mitigate 
the  energy  penalty  for  this  bond  rotation.  Clearly,  the  extra  t/-shell  valance  electron  of  tantalum 
dramatically  helps  to  stabilize  the  fault  as  compared  to  hafnium.  Since  both  systems  have 
equivalent  coordination  environments  (bond  hybridization),  the  extra  d  valance  electron  provides 
for  some  charge  delocalization  upon  rotation  and  would  explain  the  universality  for  a  stabilized 
ISF  in  all  of  the  other  group  VB  TMMCs  shown  in  the  SI  [13]. 
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This  faulted  configuration  stacks  the  metal  atoms  over  the  metal  atoms  and  the  carbon 
atoms  over  the  carbon  atoms  across  the  fault.  In  other  words,  the  atoms  directly  in  the  fault  have 
a  trigonal  prismatic  coordination  with  respect  to  the  other  species  and  tetrahedral  coordination 
with  respect  to  their  own  atom  type.  It  has  been  previously  suggested  that  in  the  synchro-shear 
mechanism  [17,25,26],  the  tetrahedral  coordination  of  the  metal  atoms  would  be  unfavorable. 
This  would  be  correct  for  HfC  but  not  necessarily  true  for  TaC,  where  an  ISF  is  favorable 
because  of  the  nature  of  its  bonds.  Finally,  it  is  interesting  to  note  that  their  appears  to  be  some 
bonding  between  carbon  atoms  in  the  TaC  ISF  though  this  is  likely  an  artifact  to  the  sensitivity 
of  the  exact  isocharge  value  used  in  visualization.  Similar  stabilization  of  ISF  energies  has  been 
noted  with  metallic  alloys  where  the  solute  atoms  contribute  excess  valence  electrons  to  the 
system  [27]. 

In  summary,  a  series  of  room  temperature  indents  in  HfC  and  TaC  were  characterized 
using  TEM  dynamical  diffraction  analysis  to  determine  the  operating  slip  planes  in  these 

materials.  We  confirmed  that  a / 2^1 10^  {110}  slip  occurs  in  HfC  and  a / 2^1 10)  {1 1 1 }  slip  in 

TaC  at  room  temperature,  with  two  orders  of  magnitude  higher  density  of  dislocations  observed 
in  TaC  than  HfC.  Using  DFT,  the  GSF  curves  for  { 1 10}  and  {111}  slip  were  computed  revealing 
that  perfect  {110}  slip  was  more  favorable  with  TaC  having  a  lower  energy  for  slip  than  HfC. 
This  result  confirmed  the  dominant  slip  system  observed  in  HfC  and  the  higher  dislocation 
density  seen  in  TaC  for  equivalent  room  temperature  indents.  The  presence  of  an  ISF  in  TaC 
promotes  the  dissociation  of  perfect  dislocations  into  partials  on  the  {111}  plane,  which  allowed 
it  to  bypass  {110}  slip.  The  extra  electron  in  TaC  as  compared  to  HfC  provides  a  more  metallic 
nature  to  the  bonds  [6,7]  enhancing  slip  and  stabilizing  the  ISF.  Coupling  prior  experimental 
work  in  the  other  group  IVB  and  VB  TMMCs  where  different  dominate  slips  planes  have  been 
reported  [1 1,14-16]  with  their  similar  GSF  curves  (see  SI),  we  have  concluded  that  the  slip 
variation  is  likely  contributed  to  the  ability  to  or  not  to  form  an  ISF.  The  stability  of  this  ISF 
appears  to  be  related  to  the  excess  d-shell  electron  between  these  two  groups  of  TMMCs.  The 
operation  of  {111}  slip  allows  for  a  larger  number  independent  slip  systems  (i.e.  5)  compared  to 
{110}  slip  (i.e.  2),  further  enhancing  ductility  of  the  material  and  changing  its  macroscopic 
properties.  This  suggest  that  engineering  the  stacking  fault  energy  to  access  or  limit  specific 
deformation  modes  through  solute  alloying  could  allow  for  the  future  design  of  tunable  hardness 
in  these  carbides  beyond  using  simple  elastic  constant  criteria  [28-30]. 
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Supplementary  Material  for:  Bonding  effect  on  the  slip  differences  in  the  B1  monocarbides 

I.  EXPERIMENTAL  DETAILS 

A.  Sample  Fabrication 

The  HfC  and  TaC  test  bars  were  produced  by  Exothennics  Inc.  using  stoichiometric 
monocarbide  powders  which  was  mechanically  pressed  and  evacuated  in  a  tantalum  container 
termed  a  'can'.  The  can  was  subjected  to  Hot  Isostatic  Pressing  (HIP)  at  200  MPa  and  1600  °C  for 
100  minutes  in  an  argon  atmosphere.  Stoichiometry  was  maintained  during  the  HIP  process,  which 
was  confirmed  (HfCo.992  and  TaCo.994)  by  LECO  combustion  and  infrared  analysis  [1]  of  the 
samples’  carbon  content.  The  post-HIP  billets  had  a  density  of  —98%  measured  by  a  water 
emersion  displacement  technique  and  compared  to  the  theoretical  densities  of  12.76  g/cc  for  HfC 
and  14.8  g/cc  for  TaC  [2]. 

B.  Indentation  and  Microscopy  Preparation 

TEM  foils  were  prepared  by  ultrasonic  ally  drilling  3  mm  discs  from  slices  of  the  HIP’ed 
billets  using  a  Fischione  Instruments  Model  170  drill  with  a  SiC  abrasive  slurry.  The  discs  were 
then  mechanically  ground  to  <  100  pm  using  SiC  and  then  indented  to  induce  room  temperature 
deformation.  A  10  x  10  matrix  of  100  Vicker’s  microhardness  indentations  using  a  300  g  load  for 
5  seconds  were  made  in  the  center  of  each  specimen  which  translates  to  a  strain  rate  of 
approximately  10"3  s"1.  This  process  of  indentation  is  similar  to  the  method  employed  by  Kim  et. 
al  [3].  In  general,  ceramics  have  been  shown  to  be  strain  rate  insensitive  for  rates  <  10+3  s"1  [4], 
For  this  test,  we  were  in  the  standard  quasi-static  strain  rate  and  our  results  should  be  universal 
over  similar  ranges  of  conditions. 
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The  un-indented  sides  were  then  dimpled  to  <  15  pm  with  6  pm  diamond  paste  polish  under 
a  rotating  15  mm  diameter  steel  wheel  using  a  Fischione  Instruments  Model  200  unit.  The  dimpled 
discs  were  Argon  ion  milled  to  electron  transparency  using  a  Gatan  Precision  Ion  Polishing  System 
(PIPS®).  The  PIPS®  dual  guns  were  operated  at  3  keV  per  gun  and  aligned  at  contrary  incident 
milling  angles  of  8°  to  achieve  perforation.  The  guns  were  operated  at  2  keV  and  6°  for  a  final 
‘clean-up’  step.  The  specimens  were  continuously  rotated  during  the  ion  milling. 

C.  Electron  Microscopy 

Conventional  dynamical  TEM  was  perfonned  with  a  FEI  F20  Tecnai  (S)TEM  operated  at 
200  keV.  The  Burgers  vectors  were  detennined  by  conventional  g  •  b  x  u  analysis,  where  g  is  the 
diffraction  vector,  b  is  the  Burgers  vector  of  the  dislocation  and  u  is  the  line  sense  of  the  dislocation 
[5].  In  Fig.  Sl(a),  the  two  groups  of  parallel  dislocations  in  the  indented  HfC  specimen  are  clearly 
visible  using  the  (200)  reflection  condition  taken  off  the  [001]  zone  axis.  Upon  tilting  to  the 
(020),  Fig.  Sl(c),  and  (131),  Fig.  S  1(f),  reflections,  these  groups  of  dislocation  become  invisible 
while  a  perpendicular  group  revealed  itself.  The  inset  images  in  Fig.  SI  are  the  two-beam 
conditions  for  each  respective  image.  The  consistent  Burgers  vector  identified  under  these 
conditions  for  the  dislocation  is  the  common  b  =  a/2(110)  type.  Considering  this  Burgers  vector 
and  other  visibility  conditions  taken  from  two-beam  conditions  off  the  [112]  and  [Oil]  zone  axes, 
the  consistent  and  stereographically  confirmed  slip  plane  was  identified  to  be  {110}  type.  In  Fig. 
Sl(h)  many  dislocations  in  the  indented  TaC  specimen  are  clearly  visible  using  the  (111) 
reflection  condition  taken  from  the  [112]  zone  axis.  The  dislocations  of  interest  are  the  long 
segments  that  ran  parallel  to  the  edge  of  the  sample.  Upon  tilting  to  the  (111)  reflection,  Fig.  S  1(1), 
the  dislocations  become  invisible.  The  consistent  Burgers  vector  identified  under  these  conditions 
for  the  dislocations  is  the  common  b  =  a/2(110)  type.  Considering  this  Burgers  vector,  other 
visibility  conditions  taken  from  two-beam  conditions  from  the  [111]  and  [Oil]  zone  axes.  The 
consistent  and  stereographically  confirmed  slip  plane  was  identified  to  be  {111}  type. 


FIG.  SI.  BF-TEM  images  with  associated  two-beam  conditions  in  the  inset.  Room  temperature 
indented  HfC  images  are  (a)-(f),  which  correspond  to  the  visibility  conditions  found  in  Table  SI. 
Room  temperature  indented  TaC  are  (g)-(l),  which  correspond  to  the  visibility  conditions  found  in 
Table  S2. 
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TABLE  SI.  “Visibility  Table”  for  HfC  tilt  experiment. 


g  (hkl) 

Possible  Burgers  vectors  (±,  times  a/2) 

Visibility 

[110] 

[110] 

[101] 

[ioT] 

[Oil] 

[Oil] 

200 

1 

1 

1 

1 

0 

0 

V 

220 

2 

0 

1 

i 

1 

1 

V 

020 

1 

1 

0 

0 

1 

1 

1 

220 

0 

2 

1 

i 

1 

1 

V 

111 

1 

0 

1 

0 

1 

0 

V 

131 

1 

2 

0 

1 

2 

1 

1 

TABLE  S2.  “Visibility  Table”  for  TaC  tilt  experiment. 
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Possible  Burgers  vectors  (±,  times  a/2) 

Visibility 

[110] 

[no] 

[101] 

[101] 

[Oil] 

[Oil] 

220 

0 

2 

1 

1 

1 

1 

V 

111 

1 

0 

0 

1 

0 

1 

V 

311 

1 

2 

1 

2 

1 

0 

V 

022 

1 

1 

1 

1 

0 

2 

V 

202 

1 

1 

0 

2 

1 

1 

V 

111 

0 

1 

1 

0 

0 

1 

1 

The  dislocation  densities,  p,  were  estimated  to  compare  of  regions  near  an  indent  versus 
regions  farther  from  an  indent.  This  was  perfonned  using  the  line  intercept  calculation  as 
outlined  by  Ham  [6]  and  by  Castro-Fernaandez  et  al.  [7]  where  n  is  the  number  of  intersections, 
L  is  the  length  of  the  grid,  and  t  is  the  thickness,  seen  in  equation  1 . 


It  would  be  reasonable  to  expect  the  dislocation  density  to  decrease  in  areas  away  from  an  indent 
for  a  brittle  material  at  room  temperature.  In  areas  near  indentations,  pmc  =  8.83  x  1013  m"2  and 
PTac  =  1 . 10  x  1014  m'2,  values  one  would  expect  for  a  worked  material.  On  the  other  hand,  in 
areas  farther  away  from  the  indents,  pHfc  =  3.12  xlO11  m"2  and  /Uac  =  1.32  xlO14  m"2,  showing  a 
marked  drop  in  dislocation  activity  for  only  HfC,  consistent  with  the  report  of  HfC  being  more 
brittle  than  TaC  under  indentation  by  Rowcliffe  et  al.  [8]. 


II.  HARD  SPHERE  MODEL 


The  hard  sphere  model  suggested  by  Van  Der  Walt  and  Sole  is  based  on  the  principle  of 
minimum  uplift.  The  amount  of  uplift,  nonnal  to  a  given  crystallographic  plane,  required  to  move 
the  atoms  in  the  Burgers  vector  direction  depends  on  the  radius  ratio  of  the  two  species  in  the 
rocksalt  structure.  The  authors  showed  that  this  would  predict  {110}  slip  for  radius  ratios  of  r/R 
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<  0.414,  {111}  slip  for  0.414<r/R<0.633  and  {001}  slip  for  r/R  >  0.633.  Note  that  the  model  also 
predicts  an  ISF  on  the  {111}  plane  when  r/R  >  0.732,  which  should  never  be  relevant  since  the 
model  would  suggest  {001}  slip.  This  model  was  correlated  to  the  slip  planes  in  ionic  rocksalt 
materials  at  later  extended  to  the  cubic  monocarbides  by  Rowcliffe  and  Hollox  [9],  however  these 
materials  are  dominated  by  covalent  and  metallic  bonding,  which  may  preclude  its  use.  Using  the 
ionic  radius  ratios  of  Toth  [10],  the  predictions  of  the  hard  sphere  model  are  shown  in  Table  S3 
along  with  the  reported  slip  planes. 


TABLE  S3.  Hard  sphere  model  slip  predictions  compared  to  observed  systems,  group  IVB  and 
VB  carbides.  Its  predictions  are  shown  to  be  incomplete.  (*Rowcliffe  [11],  +Lee  [12],  #Morgan 
[13],  &Kim  [3],  "“This  work) 


Compound 

Radius  Ratio 
r/R 

Hard  Sphere  Model 
Slip  prediction 

Observed  Slip 

r  TiC 

0.526 

111 

110*,  111* 

IVB<  ZrC 

0.483 

111 

110*+,  111+,  100+ 

^  HfC 

0.486 

111 

110*%,  m% 

f  VC 

0.576 

111 

110*,  111* 

VB  <  NbC 

0.530 

111 

110*,  1 1 1*# 

^TaC 

0.529 

111 

110%,  111*& 

III.  COMPUTATIONAL  METHODS 

A.  Generalized  Stacking  Fault  Energies 

The  simulations  of  HfC  and  TaC  were  conducted  using  DFT  as  implemented  in  the  Vienna 
ab  initio  Simulation  Package  (VASP).  A  projector  augmented  wave  (PAW)  potential  [14]  was 
used  to  describe  the  interaction  between  the  valence  electrons  and  the  ionic  cores.  The  Perdew- 
Burke-Emzerhoff  (PBE)  [15,16]  exchange  correlation  functions  within  the  generalized  gradient 
approximation  (GGA)  were  used  for  evaluation  of  the  exchange-correlation  energy.  Parametric 
testing  indicated  that  an  energy  cut-off  of  400  eV  was  sufficient  for  all  calculations  and  the  energy 
convergence  threshold  was  set  to  le-4  eV.  The  calculated  lattice  parameters  that  were  arrived  at 
using  a  Murnaghan  fit  for  the  relaxed  HfC  (4.65  A)  and  TaC  (4.47  A)  crystals  were  in  good 
agreement  with  the  experimental  values  determined  via  x-ray  diffraction  (4.64  A  and  4.47  A, 
respectively)  [17,18]  when  a  sufficiently  large  Monkhorst-Pack  k-point  scheme  was  used 
(llxllxll)  along  with  ion  relaxation  for  a  single  unit  cell. 

Using  the  lattice  parameters  determined,  supercells  with  24  to  30  atomic  layers  were 
constructed  where  the  crystallographic  plane  of  interest  was  aligned  parallel  to  the  x-y  simulation 
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axes,  (111)  and  (110),  respectively.  A  disregistry  plane  was  created  in  the  middle  of  the  supercell, 
where  all  the  atoms  above  were  rigidly  shifted.  The  GSF  energy  was  obtained  by  calculating  the 
increased  amount  of  energy  in  the  system  relative  to  the  base,  un-shifted  structure,  while  relaxing 
all  but  the  top  and  bottom  two  layers  of  atoms  in  the  supercell  with  respect  z  axis.  The  top  and 
bottom  two  layers  remained  fixed  in  all  directions.  The  k-point  meshes  were  chosen  to  be  inversely 
proportional  to  the  real  space  length  in  that  direction  such  the  ratio  of  1 1  k-points  to  the  lattice 
parameter  was  approximately  maintained.  Since  the  simulations  use  periodic  boundary  conditions, 
a  1 5  A  gap  was  placed  parallel  to  the  disregistry  plane  so  that  only  one  stacking  fault  would  be 
present. 

The  DFT  GSF  curves  for  HfC  and  TaC  are  shown  in  the  manuscript  in  Fig.  2.  The  unstable 
stacking  fault  energy  (USF)  for  the  <11 0> {110}  slip  is  2.4  and  1.7  J/m2  for  HfC  and  TaC 
respectively.  This  increases  to  4.5  and  3.4  J/m2  for  <1 1 0> { 1 1 1}  slip  in  HfC  and  TaC.  In  TaC,  the 
intrinsic  stacking  fault  energy  is  found  to  be  0.6  J/m2  and  the  USF  between  the  equilibrium  and 
ISF  was  1 .6  J/m2,  the  lowest  of  any  surface.  There  is  no  ISF  in  the  case  of  HfC  and  the  USF  is  2.8 
J/m2.  The  remaining  GSF  curves  for  the  group  IVB  and  VB  monocarbides  are  shown  in  Fig.  S2. 
The  equivalent  trends  of  an  ISF  present  in  the  group  VB  monocarbides  that  is  not  present  in  the 
group  IVB  monocarbides  is  observed  indicating  the  universal  nature  associated  with  the 
differences  in  slip  between  these  two  equivalent  symmetry  carbides. 


Group  IVB  TMMCs  Group  VB  TMMCs 


Displacement  Fraction  Displacement  Fraction 


FIG.  S2.  GSF  curves  for  the  various  group  IVB  and  VB  transition  metal  monocarbides  (TMMCs). 
Note  the  ISF  present  in  the  group  VB  series  that  is  absent  in  the  group  IVB  series. 
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B.  Peierls-Nabarro  Model 


To  model  the  dissociation  of  the  dislocation  on  the  {111}  plane,  a  standard  Peierls  Nabarro 
model  [19]  was  implemented.  The  model  treats  planar  dislocation  core  defonnation  by  accounting 
for  the  elastic  energy,  Eei,  and  misfit  energy,  Emsft,  where  the  core  structure  was  detennined  by 
minimizing  the  total  energy,  Etot .  This  relationship  is  shown  in  equations  (2)  and  (3), 

Etot  [«W]  =  Eel  T  fms/t  (2) 

Etot[u(x )]  =  -K  f*Rf*Rp(,x)p(x')\n\x  -  x'\dxdx'  +  f*Ry(u(x))dx  (3) 


where  u(x)  is  the  misfit  across  the  cut  plane,  R  is  the  cut-off  radius,  p(x)  is  the  disregistry  density, 
y(u(x))  is  the  misfit  potential,  and  K  is  a  constant  represent  the  material  properties  and  the 
dislocation  type.  For  a  screw  dislocation,  Ks  =  p/An  and  Ke  —  ix/(An(l  —  v))  for  an  edge 
dislocation,  where  p  is  the  shear  modulus  and  v  is  the  Poisson’s  ratio  for  the  material. 

By  choosing  the  trial  functions  of  equations  (4)  and  (5), 

ux\x)  =  ~(A1  arctan - h  A2  arctan - )  (4) 

7T  C-l  C2 

\  /  .  .  x-x3  .  x-x4\  b 

uz(x)  =  -  ^3  arctan  — - f  A4  arctan-^ — J  —  -  (5) 

where  b  is  the  Burgers  vector,  Xi  is  the  distance  from  the  dislocation  core,  cv  is  the  optimization 
parameter,  and  A,  is  a  proportional  constant,  and  following  the  relationship  of  equation  (6),  we  can 
also  derive  the  trial  functions  for  the  disregistry  densities,  shown  in  equations  (7)  and  (8). 


p(x)  = 


Pxto  =  -t 


A  i  Ci 


n  v(x-x1)2+cJ 


dux(x) 

clx 

+ 


a2c2 


PzM  =  -(- 


A  3C3 


■  + 


(x-x2)2+cl 

A4.C4 


(6) 

(7) 

(8) 


Notice  that  the  <x,  z>  notation  represents  the  direction.  Plane  <x,  z>  represents  a  {111}  plane  in 
the  crystal  and  x-direction  corresponds  to  a  <110>  direction  and  the  z-direction  corresponds  to  a 
<112>  direction. 


C.  Isocharge  Rendering 

The  isocharge  values  for  the  renderings  shown  in  Fig.  3.  of  the  manuscript  were  done 
using  a  value  of  0.04  e/Bohr3.  The  maximum  contour  corresponds  to  0.12  e/Bohr3  while  the 
minimum  is  that  used  to  draw  the  isosurface. 
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Abstract 

For  group  VB  transition  metal  carbides,  as  compared  to  group  IVB  carbides,  an  anomalous 
rise  in  hardness  occurs  for  substoichiometric  carbon  concentrations  as  compared  to  the 
stoichiometric  monocarbides.  Here  we  report  the  computationally  derived  elastic  constants, 
electronic  density  of  states,  and  activation  energies  for  carbon  vacancy  migration  as  a  function  of 
carbon  content  to  elucidate  their  effect  on  the  hardening  responses.  The  changes  in  elastic 
properties  with  respect  to  carbon  vacancy  concentration  show  similar  behaviors  of  elastic  softening 
and  decreasing  hardness  for  all  of  the  cubic  carbides.  The  consistent  trends  of  vacancy  diffusion 
energy  barriers  between  all  of  the  group  IVB  and  VB  transition  metal  carbides  also  suggests  that 
carbon  diffusion  may  not  be  a  significant  contributor  to  the  reported  hardness  differences. 
Consequently,  we  propose  that  the  anomalous  hardening  for  substoichiometric  behavior  is  a 
competition  between  elastic  constant  softening  and  a  microstructural-based  effect,  i.e.  domain 
hardening,  that  is  present  in  group  VB  carbides  but  not  in  group  IVB  carbides. 

Introduction 

The  group  IVB  and  VB  transition  metal  carbides  are  a  class  of  refractory  ceramics  that  are 
well  known  for  their  thermal  and  structural  stability.  These  materials  form  compounds  that  have 
the  highest  melting  temperatures  recorded  and  exhibit  exceptional  hardness,  which  are  direct 
results  of  the  mixed  covalent-metallic-ionic  bonding  present  in  these  compounds.  Thus,  they  are 
well  suited  for  structural  applications  at  high  temperatures  where  a  variety  of  strengths  are 
required. 

At  near  equal  parts  transition  metal  and  carbon,  these  materials  typically  form  the  rock  salt 
or  B 1  structure.  This  structure  is  stable  for  a  range  of  stoichiometry  allowing  a  large  amount  of 
carbon  loss.  For  the  group  VB  compounds,  the  FCC  structure  of  the  metal  atoms  is  maintained 
with  the  loss  up  to  approximately  25%  carbon  at  elevated  temperatures.  At  lower  temperatures, 
the  carbon  vacancies  can  order  resulting  the  precipitation  of  a  MeeCs  phase,  where  Me  is  the  metal 
atom  1.  Further  carbon  loss  may  result  in  the  precipitation  of  various  faulted  metal-rich  carbide 
phases  such  as  Me4C3  and  MeiC,  which  do  not  share  the  FCC  metal  sublattice  ordering.  Of  the 
group  VB  metal  carbides,  vanadium  carbide  is  unique.  It  does  not  fonn  a  stoichiometric  B1 
structure  like  TaC  or  NbC,  rather  it  phase  separates  to  substoichiometric  VCX  and  graphite,  at 
carbon  concentrations  equal  to  or  greater  than  X-0.882 .  In  contrast,  the  group  IVB  compounds 
can  sustain  much  larger  losses  of  carbon  concentrations  without  the  loss  of  the  FCC  metal 
sublattice.  Thermodynamic  and  first  principle  calculations  of  these  compounds  have  also  shown 
the  stability  of  equivalent  vacancy  ordered  phases  retaining  the  FCC  metal  lattice3,  but,  to  the 
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authors'  knowledge,  these  phases  have  not  been  experimentally  reported.  None  of  the  faulted 
metal-rich  carbides  are  known  to  fonn  in  the  group  IVB  carbides. 

It  is  well  known  that  vacancies  can  directly  affect  the  thermal  and  structural  stability  and 
exert  influence  on  mechanical  properties  of  materials.  Notably,  the  melting  temperature  of  the 
group  VB  carbides  has  a  maximum  that  is  not  located  at  stoichiometric  MeC,  but  at  some 
measurable  carbon  loss1 .  Additionally,  the  hardness  measurements  have  also  been  reported  to 
similarly  trend  with  the  peak  hardness  being  sub-stoichiometric,  Figure  1 .  For  NbCx  and  TaCxthis 
maximum  is  at  C/Me=x<0.854"7 .  While  various  studies  may  not  all  show  the  same  amount  of 
hardening  or  the  exact  same  carbon  vacancy  transitions,  it  is  clear  that  such  a  substoichiometric 
transition  occurs.  In  contrast,  studies  in  the  group  IVB  carbides  have  not  shown  a  similar 
anomalous  hardening  and  melting  temperature  behavior;  in  contrast,  in  each  of  these  properties, 
the  maximum  property  is  at  the  stoichiometric  C:  Me  ratio  of  equal  parts  7 . 


C/Me 


Figure  1:  Representation  of  hardness  as  a  function  of  transition  metal  carbide  and  carbon  content. 
The  image  is  a  rendering  of  that  reported  in  reference4 . 

Numerous  theories  have  been  developed  to  explain  this  difference  in  hardness  between  the 
group  IVB  and  VB  transition  metal  carbides,  including  optimal  filling  of  the  d-bands4 ,  changes  in 
slip  systems8  10 ,  as  well  as  the  potential  hardening  cause  by  substoichiometric  phases11 .  The  study 
by  Jhi  et  al 12  supports,  with  electronic  structure  density  functional  theory  (DFT)  calculations,  the 
notion  that  bonding  has  a  direct  effect  on  the  hardness  through  the  elastic  constants.  They  showed 
that  the  shear  modulus  of  NbCx  does  not  decrease  nearly  as  much  as  TiNx  (which  shows  a 
comparable  loss  of  hardness  similar  to  TiCx). Vacancies  have  also  been  postulated  to  control 
strength  and  plastic  flow  in  these  materials.  These  theories  7’ 13-15  assume  that  carbon  vacancy 
migration  through  the  dislocation  core  is  necessary  for  plastic  flow  at  low  and  moderate 
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temperatures.  Collectively,  these  prior  results  strongly  suggest  that  the  formation  and  migration 
of  carbon  vacancies  in  these  materials  is  critical  in  understanding  their  properties  and  ultimately 
how  to  design  these  materials. 

In  this  work,  we  investigate  the  properties  of  vacancies  in  all  of  the  group  IVB  and  VB 
transition  metal  carbides  using  electronic  structure  density  functional  pseudo-potential  theory.  We 
solely  focus  on  the  B 1  structure  with  vacancies  and  compare  our  findings  with  available  literatures. 
Notably,  we  compute  the  change  in  elastic  properties  with  respect  to  carbon  vacancy  concentration 
for  all  of  the  cubic  carbides.  The  comparison  of  these  properties  enables  us  to  determine  its 
contribution  to  the  hardening/softening  in  indentation  experiments.  We  also  compute  the  vacancy 
fonnation  energies  as  a  function  of  the  super-cell  size,  which  allows  us  to  compare  the  fonnation 
energy  and  stability  across  the  cubic  carbides.  The  computed  vacancy  migration  energy  in  all  of 
these  materials  will  be  compared  against  available  experiments  as  well  as  predict  some 
unmeasured  properties.  Finally,  we  examine  how  the  densities  of  states  in  these  materials  change 
as  a  function  of  carbon  content  in  the  B 1  structure.  While  numerous  papers  have  computed  elastic 
properties  of  the  pure  B 1  structures16"20,  none  have  comprehensively  computed  the  properties  for 
a  wide  range  of  vacancy  concentrations  or  materials.  This  self-contained  approach  will  provide 
unique  insight  into  the  similarities  and  differences  between  the  group  IVB  and  VB  cubic 
monocarbides  and  help  to  elucidate  the  anomalous  hardening  in  group  VB  carbides  that  is  not 
observed  in  group  IVB  carbides  as  shown  in  figure  1 . 

Methodology 

To  compute  the  energies  used  in  estimating  the  vacancy  fonnation  energies,  migration 
energies  and  elastic  constants,  we  employed  the  electronic  structure  density  functional  theory  as 
implemented  in  the  Vienna  Ab-Initio  Simulation  Package  (VASP)21, 22 .  For  all  of  our  calculations, 
we  use  the  Projector  Augmented  Wave  (PAW)  pseudo-potential23' 24  and  evaluate  the  exchange 
conelation  energy  using  the  Purdew-Burke-Emzerhoff  (PBE)  formulation  within  the  generalized 
gradient  approximation  (GGA)25 .  The  energy  cutoff  used  in  all  our  calculations  was  400  eV  and 
integration  in  k-space  was  performed  using  Monkhorst-Pack  with  an  integration  that  varied  with 
simulation  cell  size26. 

In  order  to  compute  properties  as  a  function  of  vacancy  concentration,  many  different 
super-cell  sizes  were  used.  The  primitive  B 1  unit  cell  consists  of  two  atoms,  but  cannot  support 
any  vacancies  without  forming  a  pure  substance.  To  resolve  this  issue,  we  used  two  types  of  super¬ 
cells:  one  that  is  based  on  an  expanded  fonn  of  the  primitive  cell  and  one  that  expanded  upon  the 
standard  FCC  unit  cell.  The  primitive  unit  cell  has  two  atoms  and  the  expanded  fonns  (equal  in 
all  directions)  have  2,  16,  54,  128,  ...  2*N3  atoms,  where  N  is  the  number  of  atoms  in  the  super¬ 
cell.  The  expanded  FCC  unit  cells  have  8,  64,  216,  512, ...  (2  N)3  atoms.  The  16  and  64  atom  unit 
cells  are  shown  in  Figure  2  for  reference.  The  advantage  of  using  both  of  these  types  of  super¬ 
cells  is  that  removing  a  single  vacancy  (either  carbon  or  metal)  results  in  the  super-cell  retaining 
its  cubic  symmetry,  even  after  relaxation,  which  allows  us  to  simulate  the  formation  of  a  relatively 
isolated  vacancy  as  well  as  determine  unambiguously  the  elastic  constants  for  the  sub- 
stoichiometric  monocarbides. 
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Figure  2:  Schematic  of  the  different  unit  cells  of  the  cubic  monocarbide  to  accommodate  the 
carbon  vacancies  without  the  loss  of  the  cubic  symmetry  (a)  Perfect  and  (b)  with  a  vacancy  for  the 
16  atoms  super-cell;  (c)  Perfect  and  (d)  with  a  vacancy  for  64  atoms  super-cell. 

Results  and  Discussion 

Elastic  constants 

Most  mechanical  properties,  i.e.  strength  and  hardness,  scale  with  elastic  constants.  We 
have  computed  both  stoichiometric,  Table  1,  and  substoichiometric,  Table  2,  elastic  properties  for 
each  of  the  carbides.  In  addition,  we  have  provided  the  lattice  constant  and  cohesive  energy  for 
each  of  the  stoichiometric  B1  structures  in  Table  1.  To  evaluate  the  enthalpy  of  formation,  the 
cohesive  energies  of  the  pure  elements,  metal  and  carbon  atoms,  were  subtracted  from  the  cohesive 
energy  of  the  compound,  with  cohesive  energies  and  lattice  constants  of  the  pure  elements  given 
in  Appendix  A  for  reference.  We  note  that  all  of  the  compounds  show  positive  enthalpies  of 
formation  and  positive  elastic  constants 


Table  1.  Properties  of  the  stoichiometric  B1  monocarbides. 


Property 

TiC 

ZrC 

HfC 

VC 

NbC 

TaC 

ao(A) 

4.331 

4.718 

4.642 

4.150 

4.501 

4.471 

Ecoh 

14.71 

15.83 

16.26 

14.01 

15.71 

17.23 

(eV/MeC) 

A  Hf  (eV) 

1.54 

1.65 

1.92 

0.787 

0.837 

1.18 

B  (GPa) 

270 

238 

253 

332 

302 

340 

Cll  (GPa) 

546 

490 

540 

669 

644 

737 

C 12  (GPa) 

133 

112 

112 

152 

129 

141 

C44  (GPa) 

162 

147 

171 

186 

166 

175 
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Table  2.The  elastic  constants  of  the  B1  mono-carbides  (in  GPa). 


Elastic 

Constant 

TiC 

ZrC 

HfC 

VC 

NbC 

TaC 

X=0.75 

Cn 

470 

403 

464 

557 

512 

579 

Cl2 

113 

105 

101 

164 

150 

172 

C44 

118 

105 

131 

107 

109 

116 

B 

233 

204 

223 

285 

270 

308 

G 

139 

120 

149 

139 

134 

147 

E 

348 

301 

366 

359 

345 

380 

X=0.875 

Cn 

520 

460 

508 

640 

608 

684 

Cl2 

117 

104 

107 

138 

126 

147 

C44 

154 

139 

162 

170 

153 

157 

B 

252 

223 

223 

305 

287 

337 

G 

171 

153 

176 

199 

181 

195 

E 

419 

374 

425 

490 

452 

488 

X=0. 96875 
Cn 

541 

484 

535 

663 

628 

709 

Cl2 

127 

108 

110 

143 

133 

147 

C44 

162 

146 

170 

184 

162 

169 

B 

265 

234 

252 

317 

287 

337 

G 

179 

161 

184 

211 

192 

207 

E 

438 

394 

445 

518 

474 

515 

To  compute  the  elastic  constants  of  the  substoichiometric  MeCx  compounds,  three  different 
super-cell  sizes:  8,  16  and  64  atoms  were  used.  The  elastic  constants  for  each  of  these  super-cells 
had  a  single  carbon  vacancy,  which  led  to  values  of  x=0.75,  0.875  and  0.96875.  The  elastic 
constants  of  the  8  and  16  atom  super-cells  without  vacancies  were  also  computed  for  comparison 
and  showed  negligible  variation  from  the  values  in  Table  1.  Though  the  extreme  carbon  loss  can 
result  in  the  precipitation  or  co-precipitation  of  metal-rich  carbide  phases,  these  simulation  results 
provide  the  ability  to  clearly  probe  the  trend  in  the  B1  compounds’  elastic  constants  with  carbon 
loss,  something  that  can  be  difficult  to  do  experimentally.  The  cubic  elastic  constants  are  listed  in 
Table  2  for  all  of  the  substoichiometric  compounds. 

The  Hill  average  elastic  constants,  which  are  an  average  of  the  Voigt  and  Reuss  elastic 
constants  27 ,  given  below,  were  detennined  from  the  computed  cubic  constants.  This  methodology 
allows  the  trends  in  elastic  constants  to  be  easily  compared.  The  Voigt  and  Reuss  shear  modulus 
are  defined  by28 

Gv  =  i(C11-C12  +  3C44)  (1) 

Gr  =  5[4(S11-S12)  +  3S44]1  (2) 

And  the  Hill  average  shear  modulus  defined  as 

G  =  i(Gv  +  GR)  (3) 
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The  Voigt  and  Reuss  bulk  modulus  are  the  same  as  the  bulk  modulus  of  the  cubic  crystal,  and 
hence  so  is  the  Hill  bulk  modulus.  The  other  elastic  constants  can  then  be  easily  derived  from  the 
shear  and  bulk  modulus. 

The  average  shear  modulus  for  all  six  cubic  carbides  are  plotted  in  Figure  3  as  a  function 
of  the  carbon  content.  Figure  3(a)  shows  a  direct  comparison  of  the  absolute  average  shear 
modulus  value  of  all  six  compounds.  Notably,  the  group  VB  carbides  are  stiffer  than  the  group 
IVB  carbides.  Figure  3(b)  is  the  same  data  normalized  by  the  shear  modulus  at  zero  vacancy 
concentration.  From  this  plot,  it  is  obvious  that  all  of  the  carbides  appear  to  elastically  soften  with 
carbon  vacancy  concentration  at  the  same  rate.  The  only  significant  spread  occurs  at  x=0.75,  which 
may  be  a  result  of  the  high  constraints  on  the  atoms  in  these  highly  symmetric  (and  small)  super 
cells. 


(a)  (b) 


X=C/Me  X=C/Me 

Figure  3:  (a)  Computed  average  shear  modulus  as  a  function  of  carbon  content  for  the  various 
transition  metal  carbides  (b)The  same  shear  modulus  normalized  by  the  shear  modulus  at  x  = 
1.0. 


The  hardness  for  all  the  carbides  was  estimated  using  these  computed  elastic  constants 
according  to  29 

Hv  =  2(k2G)0585-3  (4) 

where  k  is  Pugh's  modulus,  which  is  the  ratio  of  the  shear  modulus,  G,  and  bulk  modulus,  B.  These 
results  are  plotted  in  Figure  4.  In  general,  the  hardness  increases  with  carbon  content.  The 
experimentally  reported  anomalous  hardness  peaks  for  substoichiometric  TaCx  and  NbCx  near 
X~0.844was  not  captured  in  these  simulations.  This  demonstrates  that  the  anomalous  hardness  is 
not  a  direct  consequence  of  the  trends  associated  with  the  elastic  constants. 
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X=C/Me  X=C/Me 

Figure  4:  (a)  The  predicted  hardness  for  each  of  the  cubic  carbides  and  (b)  the  same  hardness 
nonnalized  by  the  hardness  at  stoichiometry,  i.e.  x  =  1.0. 

Electronic  Structure  and  Bonding 

Some  works  have  suggested  that  the  bond  structure  of  the  carbides  contributes  directly  to 
the  anomalous  hardness  observed  in  TaC  and  NbC4' 1 .  The  elastic  constants  computed  above 
provide  direct  evidence  to  the  contrary;  all  the  transition  metal  carbides  elastically  soften  at 
approximately  the  same  rate  with  carbon  loss.  The  electronic  density  of  states  (DOS)  in  the  MeC, 
MeCo.875  and  MeCo.75  structures  for  all  of  the  cubic  transition  metal  carbides  has  been  computed. 
Examples  of  the  total  DOS  computed  for  TaC  (group  VB  carbide)  and  TiC  (group  IVB  carbide) 
are  shown  in  Figure  5.  The  DOS  peak  for  TaC  is  near  -5  eV,  Figure  5(a),  and  is  associated  with 
the pd  bonding  between  the  tantalum  and  carbon  atoms,  which  is  covalent  in  nature.  It  appears  that 
this  peak  decreases  with  vacancy  concentration,  signaling  a  decrease  in  the  contribution  from  these 
bonds.  Similarly,  a  slight  increase  in  the  DOS  at  the  Fenni  level  accompanies  this  decrease  which 
is  interpreted  as  an  increase  in  the  dd  bonding  between  the  tantalum  atoms.  The  TiC  pd  and  dd 
bonding  shows  a  similar  trend,  Figure  5(c),  although  the  details,  such  as  the  peak  location  of  the 
pd  bonding  ( — 3  eV)  is  slightly  different  than  that  in  TaC,  Figure  5(a).  These  trends  in  the  DOS 
with  respect  to  carbon  content  are  consistent  between  all  of  the  other  transition  metal  carbides. 
This  transition  from  covalent  towards  metallic  character  results  in  a  decrease  in  the  strength  of  the 
bonds  which  is  manifested  in  the  elastic  constants  reported  in  Tables  1  and  2. 
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(a) 


(b) 


Figure  5:  (a)  The  total  density  of  states  for  TaC,  TaCo.875  and  TaCo.75.  The  DOS  is  normalized  by 
the  number  of  formula  units.  Alternatively,  the  normalization  can  be  thought  of  as  the  number  of 
metal  atoms,  (b)  A  comparative  DOS  plot  of  TiC. 

Activation  energies  for  diffusion. 

Recognizing  that  carbon  vacancy  diffusion  may  contribute  to  plastic  flow14,  30  32  ,  the 
anomalous  peak  in  hardness  may  be  associated  with  carbon  diffusion.  It  has  been  proposed  that 
carbon  vacancy  diffusion  is  the  rate  limiting  step  in  plastic  deformation.  This  is  because  the 
activation  energy  associated  with  plastic  deformation  is  reported  to  be  close  to  that  of  carbon 
vacancy  diffusion33  .For  example,  at  low  temperatures,  TaC  is  reported  to  slip  on  {111}  planes  9- 
34  and  is  this  is  same  slip  planes  for  most  cubic  carbides  at  elevated  temperatures  8 .  Slip  on  these 
close-packed  planes  is  assisted  by  the  dissociation  of  the  perfect  dislocation  into  partial 
dislocations14, 35, 36.  This  splitting  motion  can  be  aided  if  the  carbon  atoms  diffuse  through  the 
dislocation  core,  or  shear  in  a  related  fashion,  such  that  metal  atoms  do  not  glide  over  other  metal 
atoms,  a  plastic  deformation  mechanism  termed  the  synchro-shear  mechanism  1S.  It  is  generally 
accepted  that  carbon  diffusion  is  of  particular  importance  to  the  deformation  mechanics  of  these 
carbides. 

Consequently,  simulations  were  conducted  on  the  formation  energy  and  migration  energy 
of  carbon  atoms  for  both  carbide  groups.  Since  carbon  vacancies  readily  form  in  the  transition 
metal  carbides1,  this  implies  that  the  energy  to  form  a  carbon  vacancy  must  be  relatively  low. 
Several  studies  using  electronic  structure  theory  have  shown  that  the  vacancy  formation  energies 
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are  indeed  low  in  TiC  and  NbC37.  However,  many  of  these  calculations  are  done  in  relatively 
small  super-cell  volumes  and  wider  comparisons  between  other  group  IVB  and  VB  carbides  have 
not  been  made.  Here,  we  computed  the  energy  to  remove  a  carbon  vacancy  from  a  number  of 
super-cell  sizes  that  range  between  4  and  108  metal  atoms.  The  energy  of  the  super-cell 
configuration  with  a  single  carbon  vacancy  is  compared  to  a  super-cell  without  a  carbon  vacancy 
to  determine  the  energy  to  remove  the  carbon  atom.  To  compute  the  vacancy  formation  energy, 
one  must  know  the  chemical  potential  of  a  carbon  atom  in  the  carbide.  Given  that  a  definitive 
measure  of  the  chemical  potential  does  not  exist  in  this  state,  we  have  used  the  cohesive  energy  of 
a  carbon  atom  in  graphite.  While  this  choice  may  not  be  ideal,  it  has  provided  substantial  insight. 
If  we  define  the  vacancy  fonnation  energy  in  tenns  of  the  number  of  metal  atoms,  N,  in  the  super¬ 
cell,  the  vacancy  fonnation  energy  is  given  as: 

Ev(MeC)  =  E(NMe,  NC)-E(NMe,  (N-l)C)-E(graphite)  (5) 

and  the  super-cell  size  was  increased  until  the  energy  converged.  This  converged  value  is  an 
approximation  to  the  vacancy  formation  energy  in  the  stoichiometric  monocarbide.  As  is  pointed 
out  by  Tan  et  al. 37 ,  we  can  also  define  an  average  vacancy  fonnation  energy  in  MeCx  as 
E*(MeCx=(N_1)/N)  =  E(NMe,  NC)-E(NMe,  (N-l)C)-E(graphite)(6) 

The  vacancy  fonnation  energy  for  a  metal  atom  is  defined  in  a  similar  manner.  However,  in  this 
case  the  chemical  potential  of  the  metal  atom  is  taken  as  the  cohesive  energy  of  the  pure  metal. 

The  vacancy  fonnation  energy  for  both  carbon  and  metal  vacancies  converge  rapidly  with 
cell  size,  as  shown  in  the  representative  plot  for  hafnium  carbide  in  Figure  6(a).  This  plot  not  only 
demonstrates  the  rapid  convergence  of  the  vacancy  formation  energy  with  the  number  of  atoms  in 
the  super-cell  but  also  highlights  the  asymmetry  of  the  vacancy  formation  energies  between  the 
metal  and  carbon  vacancies.  Obviously,  it  is  much  easier  to  form  an  isolated  carbon  vacancy  than 
a  metal  vacancy.  Table  3  lists  the  converged  carbon  and  metal  vacancy  formation  energies  for 
stoichiometry  for  each  of  the  carbides. 


Table  3:  The  formation  energy  (in  eV)  of  the  carbon  and  metal  vacancies  computed  from  the 
largest  super-cell. _ 


Material 

Carbon  Vacancy 

Metal  Vacancy 

TiC 

0.42 

8.6 

ZrC 

0.28 

9.4 

HfC 

1.1 

9.3 

VC 

-0.93 

4.5 

NbC 

-0.42 

4.1 

TaC 

0.26 

3.5 
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Figure  6:  (a)The  metal  and  carbon  vacancy  formation  energy  in  hafnium  carbide  plotted  as  a 
function  of  the  total  atoms  in  the  super-cell. (b)The  average  carbon  vacancy  formation  for  all  of 
the  compounds  as  a  function  of  the  carbon  concentration.(c)The  carbon  vacancy  migration  energy 
plotted  as  a  function  of  carbon  content.  The  filled  markers  represent  results  from  calculations  with 
a  rigid  cell  while  the  open  markers  are  results  where  the  cell  is  allowed  to  change  shape  and  size. 

The  group  IVB  TiC,  ZrC,  and  HfC  as  well  as  the  group  VB  TaC  all  have  positive  carbon 
vacancy  formation  energies.  This  implies  that  the  stoichiometric  B1  phase  is  stable,  as 
experimentally  observed.  The  negative  value  for  VC  would  indicate  that  the  B1  phase  is  unstable 
or  would  spontaneously  form  vacancies.  Recall,  VC  has  not  been  reported  to  form  the 
stoichiometric  B1  phase2  and  these  results  would  agree  with  experiments.  In  fact,  most  phase 
diagrams  of  vanadium  carbide  have  shown  the  carbide  separating  between  VCX  and  graphite  near 
X=0.88.  Note  that  the  average  vacancy  formation  energy  becomes  positive  around  x=0.85,  Figure 
3(b). 


The  NbC  was  also  noted  to  have  a  negative  vacancy  formation  energy  value,  Table  3,  yet 
is  known  to  experimentally  fonn  the  stoichiometric  B 1  phase38 .  The  negative  value  arises  because 
of  the  chemical  potential  used  is  the  DFT  computed  cohesive  energy  of  graphite,  i.e.  7.85  eV/atom, 
and  as  noted  previously  would  not  be  the  chemical  potential  of  carbon  explicitly  in  these  carbides. 
Consequently,  this  has  affected  the  accuracy  of  some  of  the  values  reported.  For  the  experimental 
cohesive  energy  of  graphite,  this  difference  is  approximately  -0.45  eV39  to  the  DFT  computed 
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value;  this  discrepancy  is  commonly  known  to  occur  in  DFT  calculations  for  pure  carbon.  If  the 
experimental  value  is  used,  the  vacancy  fonnation  energy  for  NbC  becomes  essentially  zero  and 
the  B 1  phase  is  stable.  Using  the  experimental  cohesive  energy  value  for  graphite  will  shift  all  the 
computed  vacancy  formation  energies  up  by  0.45eV  in  Table  3.  Since  we  are  consistently  using 
DFT  computed  values  throughout  this  work,  we  have  chosen  to  retain  those  values  here  for 
consistency. 

Coupling  the  computed  vacancy  formation  energy  with  the  migration  energy,  the  activation 
energy  for  carbon  vacancy  diffusion  can  be  now  determined.  The  migration  energy  was  computed 
assuming  that  the  diffusion  occurs  along  the  <11 0>  directions  between  metastable  states  for  single 
carbon  vacancies  created  above.  The  minimum  energy  path  is  detennined  using  the  nudged- 
elastic-band  (NEB)  method40  to  relax  the  initial  path,  which  is  a  linear  interpolation  between  the 
initial  and  final  configuration.  We  used  eleven  discrete  images  in  our  calculation  along  this  initial 
and  final  states  with  the  migration  energy  determined  for  carbon  contents  of  x=0.75,  0.675,  0.969, 
and  0.99.  To  determine  the  effects  of  constraints,  the  migration  energy  for  fixed  cell  and  flexible 
geometries  were  employed.  For  the  fixed  geometry,  the  super-cell  is  not  allowed  to  change  shape 
or  size  during  NEB  minimization  while  both  are  allowed  to  change  for  the  flexible  geometry.  As 
one  could  expect,  this  geometry  difference  causes  variations  in  the  computed  migration  energies 
plotted  in  Figure  6(c),  with  the  smaller  cell  geometries  (or  lower  carbon  concentrations)  exhibited 
the  most  significant  deviations.  If  the  cell  geometry  is  held  fixed,  the  migration  energies  tended  to 
decrease  because  the  lower  carbon  concentrations  have  smaller  lattice  constants,  hence  bond 
distances,  making  migration  more  difficult.  Interestingly,  the  lower  migration  energy  of  the 
flexible  cells  at  lower  carbon  concentrations  suggests  that  these  cells  may  intrinsically  be  more 
flexible  for  the  diffusion  of  atoms. 

As  a  result  of  the  migration  energy  spread  in  Figure  6(c),  particularly  at  the  lower  carbon 
vacancy  concentrations,  it  is  not  immediately  obvious  which  value  should  be  taken  as  the  correct 
migration  energy.  This  computational  spread  may  provide  insight  into  the  experimental  scatter 
reported  for  the  activation  energies  in  the  various  transition  metal  carbides,  Table  4,  because  it 
may  not  have  a  fixed  value.  Since  the  computed  values  near  stoichiometry  are  close  to  the  mean 
of  the  values  at  higher  carbon  vacancy  concentrations,  these  results  do  not  support  a  concentration 
dependent  migration  energy.  Therefore,  we  have  assigned  average  migration  energy  to  be  the  value 
near  stoichiometry. 
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Table  4.The  formation  and  migration  energy  (in  eV)  of  the  carbon  vacancies  and  experimentally 
detennined  activation  energies.  The  values  in  parenthesis  in  the  experiments  are  the  carbon 
concentration. 


Material 

Ff 

cv 

E“ 

Q 

Experiments 

TiC 

0.42 

3.7 

4.1 

4.13  (0.97),  4.77 
(0.67),  2.15  (0.67), 

4.63  (0.89),  4.81 
(0.47),  4.25  (0.87- 
0.96)47'50 

ZrC 

0.28 

4.3 

4.6 

4.91  (0.92)  ,4.94 
(0.97)51 

HfC 

1.1 

4.8 

5.9 

VC 

-0.93 

3.0 

3.0 

2.91  (0.75),  3.11 
(0.92),  3.69  (0.84)51, 50 

NbC 

-0.42 

3.6 

3.6 

3.8  (0.7-0.93),  4.08- 
4.34  (0.75-0.97) 

TaC 

0.26 

4.0 

4.3 

3.93  (0.99),  3.73 

3.73,  4.3,  3. 852 53 

Summing  the  migration  and  formation  energy,  the  computed  activation  energy  for 
diffusion  is  tabulated  in  Table  4.  For  NbC  and  VC,  who’s  computed  formation  energies  are 
negative,  we  have  taken  the  formation  energy  in  each  of  these  two  special  cases  to  be  zero. 
Considering  the  results  in  Table  4,  there  is  no  obvious  difference  in  activation  energies  for 
diffusion  that  could  contribute  to  plastic  flow  between  the  group  IVB  and  VB  transition  metal 
carbides  with  respect  to  carbon  concentration  that  would  explain  the  anomalous  hardness. 


Phase  and  Microstructure 


Considering  those  elastic  constants,  direct  bonding,  and  vacancy  properties  are  not  direct 
causes  of  the  anomalous  hardening,  an  alternate  explanation  is  necessary.  The  anomalous  hardness 
is  only  reported  to  be  near  X-0.864, 7,  which  is  near  the  composition  for  the  MeeCs  phase.  Of  the 
transition  metal  carbides,  this  phase  has  only  been  experimentally  reported  only  in  the  group  VB 
carbides  where  the  anomalous  hardness  is  reported  in  the  tantalum  and  niobium  carbides.  Any 
explanation  that  relates  the  anomalous  hardness  to  the  MceCs  phase  must  account  for  its  absence 
of  a  hardness  peak  in  vanadium  carbide.  This  can  be  accounted  for  considering  that  the  phase  and 
microstructures  that  form  around  these  compositions. 

In  Figure  4,  the  computed  hardness  of  the  rock  salt  phase  with  vacancies  demonstrates  a 
universal  decrease  in  hardness  with  carbon  loss.  In  Figure  7,  we  have  plotted  the  hardness 
associated  with  specific  crystal  structures  that  are  associated  with  specific  stoichiometries,  i.e. 
MeeCs,  Me4C3,  etc.  using  prior  reported  elastic  constants  for  those  phases3'41,44,45 .  This  plot,  similar 
to  Figure  4,  reveals  that  the  theoretical  elastic  constant  derived  hardness,  equation  (4),  decreases 
with  depleted  carbon.  Clearly,  the  formation  of  a  specific  metal-rich  carbide  phase  also  cannot 
directly  explain  the  anomalous  hardness.  Rather,  we  believe  the  combination  of  phase  and 
microstructure  are  necessary  to  elucidate  this  behavior,  as  only  elastic  constant  contributions,  not 
microstructure,  are  captured  in  the  computed  hardness  values  using  equation  (4). 
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Figure  7:  The  predicted  hardness  from  the  elastic  constants  (reference  given  after  each)  for  the 
specific  phase  present  at  that  composition.  (a)Tantalum  carbides44  (b)  Niobium  carbides41 
(c) Vanadium  carbides45  and  (d)  Hafnium  carbides3 . 

Following  the  ideas  of  Williams42  and  Morgan  and  Lewis11,  the  peak  in  hardness  is  related 
to  vacancy  ordering  in  TaCx  and  NbCx42 .  Morgan  and  Lewis11  argued  that  the  rise  in  hardness  in 
NbCx  occurs  due  to  the  superposition  of  two  mechanisms:  the  hardening  from  domain  boundaries 
in  the  ordered  NbeCs  phase  and  the  increased  Peierls  stress  due  to  the  increase  in  cohesive  energy 
of  the  ordered  crystal.  In  their  work,  the  isolated  contributions  of  elastic  constants  and  vacancy 
properties  were  not  individually  addressed,  as  done  here,  so  that  a  firm  conclusion  of  these 
contributions  could  not  be  eliminated  as  a  possibility.  Domain  hardening  in  V6C5  has  been  reported 
to  follow  essentially  a  Hall-Patch  relationship  43 ,  which  we  expect  to  be  true  in  other  similar  group 
VB  tantalum  and  niobium  carbides,  though  this  experimental  work  has  not  yet  been  reported. 
Given  that  the  MeeCs  phase  is  elastically  softer  than  MeC,  Figure  7,  the  combined  microstructure 
of  ordered  domains  within  the  host  rock  salt  grains  must  be  contributing  to  the  anomalous  increase 
in  hardness.  Considering  our  prior  findings  in  section  3.1,  the  contribution  from  the  cohesion  of 
the  solid  must  be  a  weak  effect  as  the  elastic  constants  continually  decrease  with  a  decrease  in 
carbon  content.  Hence  the  dominant  influence  of  the  ordered  MceCs  phase  must  be  significant  in 
increasing  the  hardness  in  the  microstructure  via  domain  hardening. 
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The  lack  of  anomalous  hardness  in  vanadium  carbide  is  not  caused  by  a  lack  of  domain 
hardening;  contrary,  it  is  the  only  carbide  system  of  these  groups  that  has  experimentally 
demonstrated  domain  hardening  43 .  But  unlike  the  other  carbides,  this  material  is  unable  to 
thermodynamically  stabilize  the  stoichiometric  monocarbide,  thus  vanadium  carbide  has  likely 
peaked  its  hardness  at  V6C5  but  there  is  no  other  phase  above  this  concentration  to  demonstrate  the 
loss  in  hardness. 

The  lack  of  anomalous  hardness  in  the  group  IVB  carbides  is  contributed  to  the  inability 
of  these  phases  to  precipitate  out  the  vacancy  ordered  MeeCs  phase.  The  lack  of  this  phase  results 
in  no  opportunity  for  domain  hardening;  hence  the  hardness  is  solely  dictated  by  the  change  in 
elastic  constants. 

Summary  and  Perspective 

In  this  paper,  the  elastic  constants,  electronic  density  of  states,  and  the  activation  energy 
for  carbon  vacancy  migration  have  computed  for  the  group  IVB  and  VB  carbides  as  a  function  of 
carbon  content.  The  elastic  constants  showed  a  monotonic  decrease  with  carbon  loss.  This 
translated  to  a  decrease  in  covalent  to  more  metallic  type  bonding.  The  relative  constant  values  of 
the  activation  energy  for  diffusion  with  respect  to  carbon  concentration  suggest  that  carbon 
diffusion  and  its  contribution  to  plastic  flow  and  hardness  cannot  account  for  the  observed 
anomalous  hardness.  The  apparent  correlation  between  the  composition  of  the  peak  in  the 
anomalous  hardness  and  the  MeeCs  phase  strongly  supports  the  idea  that  this  phase  is  contributing 
to  the  anomalous  hardness.  This  is  strengthened  by  the  fact  that  the  anomalous  hardness  is  only 
observed  in  materials  that  experimentally  form  this  phase.  However,  our  computational  work 
demonstrates  that  the  carbides  should  all  elastically  soften  with  carbon  loss  and  that  the  anomalous 
hardness  cannot  solely  be  contributed  to  the  elastic  constants  of  this  phase.  Rather  its  precipitation 
as  domains  in  the  rock  salt  matrix  grains  provides  for  a  microstructural  hardening  effect,  i.e. 
domain  hardening (43).  The  computational  study  here  has  provided  further  insights  that  the  cohesion 
of  the  solid  is  less  significant  in  this  contribution  to  hardness  as  previously  speculated  4- n- 12.  Thus, 
the  anomalous  hardness  is  a  result  of  a  competition  between  elastic  constant  softening  and 
microstructure-based  hardening  mechanisms. 

Acknowledgments 

XX  Yu  and  GB  Thompson  recognize  Air  Force  Office  of  Scientific  Research  grant  FA9550- 
12-1-0104,  Dr.  Ali  Sayir  program  manager  and  the  University  of  Alabama  Office  for  Research 
Post-Doctoral  Fellowship  Program.  Professor  C.H.  Turner  of  the  University  of  Alabama  and  the 
Alabama  Supercomputer  (http://www.asc.edu/)  are  also  acknowledged  for  access  to  their 
computational  resources. 


Appendix 

Properties  of  the  Pure  Elements 

The  group  IVB  pure  elements  form  a  hexagonal  close-packed  structure  and  the  group  VB 
from  body-centered  cubic  structures.  The  lattice  constants  and  cohesive  energies  were  computed 
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for  these  structures  using  the  aforementioned  DFT.  The  group  VB  elements  were  computed 
previously  by  the  authors  in  a  separate  study  but  are  reproduce  here  for  convenience  46 .  The 
cohesive  energy  and  lattice  constants  for  both  graphite  and  diamond  were  computed  as  references 
for  carbon.  As  noted  in  the  main  text,  the  cohesive  energy  for  both  graphite  and  diamond  are  about 
0.4  eV  too  high.  Also,  the  c  lattice  parameter  of  graphite  is  grossly  incorrect  as  DFT  has  difficulty 
in  capturing  the  out-of-plane  sigma  bonds.  Both  these  deficiencies  are  well  known. 


Table  A.  1 :  Properties  of  the  pure  elements 


Property 

Ti 

Zr 

Hf 

V 

Nb 

Ta 

Graphite 

Diamond 

a  (A) 

2.93 

3.23 

3.19 

2.98 

3.33 

3.31 

2.47 

3.57 

c  (A) 

4.62 

5.18 

5.05 

7.72 

Ecoh  (eV) 

5.31 

6.33 

6.48 

5.37 

6.93 

8.20 

7.85 

7.72 
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Abstract 

Using  a  variable-composition  ab  initio  evolutionary  algorithm,  the  stability  of  various  tantalum 
carbide  compounds  at  ambient  pressure  and  at  OK  was  investigated.  The  results  revealed  that  TaC, 
TaeCs  and  Ta2C  are  the  lowest  energy  configurations  with  Ta4C3  and  TaXT  having  slightly  higher 
energies.  The  vacancy  ordered  TaeCsphase  had  three  energetically  degenerate  structures.  A 
competition  between  vacancy  ordered  and  stacking  fault  variation  of  the  phases  was  seen,  with  the 
latter  becoming  more  favorable  with  lower  carbon  content.  The  close  formation  enthalpy  of  each 
stable  and  metastable  phase  appears  to  ‘frustrate’  the  carbide  in  the  co-precipitation  of  multiple 
phases  for  substoichiometric  compositions.  Density  functional  theory  calculations  also  provided 
the  elastic  constants  for  each  of  the  stable  and  metastable  phases.  As  the  carbon  content  increased, 
the  elastic  constants  values  increased.  This  has  been  associated  with  the  change  in  metallic  to  more 
covalent  bonding  of  the  carbide  from  the  density  of  states.  The  collective  results  of  this 
computational  work  provide  insight  into  why  specific  tantalum  carbide  phases  form  and  the 
consequences  they  have  on  microstructure  and  properties. 

Introduction 

Tantalum  carbides,  along  with  others  in  the  group  VB  transition  metal  carbides,  exhibit  a  range 
of  high  to  ultrahigh  melting  temperatures,  microstructures  controlled  by  the  co-precipitation  of 
metal-rich  carbide  phases,  and  varied  mechanical  behaviors  dependent  on  the  phase  content 
[1,2, 3, 4, 5, 6].  Consequently,  these  materials  have  been  proposed  for  various  thennal  heat 
protection,  automotive  wear  resistant  liners,  and  other  types  of  thermo-mechanical  loading 
applications  [7,8,9].  Though  these  carbides  maybe  classified  as  ceramics,  which  are  often 
considered  to  be  hard  and  brittle  materials,  some  tantalum  carbide  phases  exhibiting  significant 
plasticity  that  exceeds  30%  at  elevated  temperatures  [3,10].  Thus,  the  potential  exists  in  this  system 
to  create  multiphase  carbides  that  exhibit  varied  mechanical  properties,  which  provides  unique 
opportunities  to  tailor  their  microstructure-mechanical  property  relationships.  For  example,  it  has 
been  reported  that  a  high  volume  fraction  of  Ta4C3  within  a  tantalum  carbide  matrix  increases  the 
room  temperature  Kic  to  ~  12- 14  MPaVm  [10,1 1],  though  the  associated  toughening  mechanisms 
for  this  increase  have  not  been  fully  elucidated  but  is  likely  linked  to  the  material’s  microstructure. 
These  high  fracture  toughness  microstructures  consist  of  several  co-precipitated  nanoscale  Ta4C3 
laths,  which  are  believed  to  contribute  to  extrinsic  toughening.  The  precipitation  and  orientation 
relationship  of  the  Ta4C3  phase  is  dictated  by  the  close  packed  planes  and  directions  of  the  parent 
material  [1,12,13].  If  the  parent  phase  is  TaC,  with  four  {111}  close-packed  plane  variants,  a 
crisscross  pattern  of  Ta4C3  laths  forms  in  an  equiaxed  grain  structure,  figure  1(a).  In  contrast, 
precipitation  of  Ta4C3  from  a  parent  Ta2C  phase  results  in  parallel  lath  fonnation  because  of  Ta2C’s 
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single  close  packed  basal  plane  [2],  figure  1(b).  These  grains  also  appear  to  be  more  acicular  [2], 
Unlike  the  group  VB  metal  carbides,  the  group  IVB  metal  carbides  do  not  form  these  metal-rich 
phases  [14]  limiting  the  opportunities  to  tailor  the  microstructure,  hence  mechanical  properties, 
through  phase  transformation  routes.  By  understanding  the  phase  equilibrium  of  the  group  VB 
metal  carbides  insights  can  be  gained  that  are  is  crucial  in  tuning  the  microstructures  and 
mechanical  properties  of  these  unique  materials. 

TaC  +  Ta4C3  TaC  +  Ta4C3  +Ta2C 


FIGURE  1:  (a)  Micrograph  of  TaCx  with  crisscross  laths  of  Ta4C3  (pointed  to  by  arrows)  for  a 
55Ta:45C  sample,  (b)  Micrograph  of  TaCx  +  Ta4C3  +  TaiC  for  a  63Ta:37C  sample.  The  x-ray 
diffraction  of  these  phases  with  these  microstructures  can  be  found  in  reference  [1]. 

In  this  paper,  we  explore  the  stability  and  structure  of  the  tantalum  carbide  phase  diagram  using 
electronic  structure  density  functional  theory  calculations.  To  date,  there  has  not  been  a 
comprehensive  and  inclusive  computational  investigation  of  the  different  phases  that  comprise  the 
tantalum  carbon  system  [15].  In  contrast,  there  has  been  some  limited  computational  investigations 
into  the  group  VB  niobium  and  vanadium  carbides  [16,17].  According  to  the  phase  diagram  of  the 
Ta-C  system  provided  by  Gusev  et  al.  [18,19],  there  are  four  compounds  in  the  tantalum  rich 
portion  of  the  diagram:  TaC,  TaeCs,  Ta4C3  and  Ta2C.  The  structure,  space  group,  and  lattice 
parameter  for  these  phases  are  tabulated  in  Table  I.  The  monocarbide  phase  of  TaC  forms  the  B1 
or  rock  salt  crystal  structure  and  similar  to  other  group  VB  metal  carbides  is  not  a  line  compound. 
It  can  accommodate  the  loss  of  some  carbon  creating  several  site-occupancy  defects  [20]  with  an 
accompanying  change  in  lattice  parameter  [18].  For  substoichiometric  TaCi-x  at  temperatures 
lower  than  approximately  1500  K,  a  TaeCs  phase  can  precipitate.  This  phase’s  experimental 
presence  is  very  uncommon  and  has  been  suggested  to  require  a  very  low  cooling  rate  (-0.15 
K/minute)  to  form  [19].  Consequently,  it  is  not  even  shown  on  the  Ta-C  phase  diagram  by  Fisk 
[21].  Gusev  etal.  [14,18]  proposed  that  this  phase  is  derived  from  cubic  TaC  i-x  where  the  vacancies 
create  a  vacancy  ordered  structure  that  is  either  monoclinic  or  trigonal.  Unlike  the  tantalum 
carbides,  this  MceCs  phase  (where  Me  is  a  metal)  is  more  frequently  observed  in  the  niobium  and 
vanadium  carbides  [18],  which  are  also  group  VB  metal  carbides.  The  phase  stability  between  the 
sub-stoichiometric  TaCi-x  structure  and  these  yet  to  be  definitively  determined  TaeCs  structures 
provides  clear  motivation  for  further  computational  investigation. 

Unlike  the  TaeCs  phase,  the  rhombohedral  C-Ta4C3  phase  has  been  experimentally  observed 
by  several  groups  [19], [22], [23, 24].  Similar  to  the  TaeCs  phase,  there  has  been  some  debate  on  its 
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thermodynamic  stability.  In  the  pioneering  work  of  Rudy  and  Hannon  [12],  they  noted  that  the 
Ta4C3  phase  formed  under  compressive  stresses  by  epitaxially  growing  off  of  the  TaCi-x  phase 
when  the  temperature  was  too  low  to  phase  transform  directly  into  TaiC.  After  400  hours  at 
1700°C,  Ta4C3  decomposed  suggesting  to  them  that  Ta4C3  phase  was  metastable.  Considering  this 
phase’s  significant  effect  on  tantalum  carbide’s  mechanical  behavior  as  noted  above  [11], 
understanding  its  thermodynamic  stability  is  essential.  This  phase  has  also  been  shown  to 
experimentally  co-exist  in  the  presence  of  TaC  and  TaiC;  yet  thermodynamics  dictates  that  only 
two  phases  should  exist  in  binary  systems  providing  yet  further  need  to  elucidate  the  phase 
equilibrium  between  all  of  these  phases. 

Further  reduction  in  carbon  will  stabilize  the  Ta2C  phase  where  the  metal  atoms  form  a 
hexagonal  closed-packed  arrangement.  This  stoichiometric  ratio  compound  exhibits  an  order  - 
disorder  transition  with  temperature.  At  temperatures  less  than  -1800  K,  the  a-Ta2C’s  carbon 
sublattice  is  ordered,  occupying  every  other  plane  of  octahedral  interstices,  while  at  elevated 
temperatures  this  sublattice  disorders  [25].  Further  loss  in  carbon  does  not  result  in  any  additional 
carbide  compounds.  The  metallic  body  centered  cubic  Ta  phase  exhibits  negligible  carbon 
solubility  except  at  temperature  near  3000  K  [19]. 

Tantalum  carbides  represent  a  rich  collection  of  phases  where  the  stability  and 
microstructure  are  still  not  completely  understood.  The  clear  links  of  phase  content  with 
microstructure  formations  and  mechanical  properties  dictates  the  need  for  a  definitive  equilibrium 
investigation.  This  work  aims  to  clarify  the  energetics  associated  within  the  tantalum  carbide 
compounds  to  better  understand  their  intrinsic  stability. 

Computational  methodology 

The  stoichiometric  compounds  and  their  crystal  structures  were  predicted  by  using  the 
evolutionary  algorithm  implemented  in  the  USPEX  code  [26,  27]. This  approach  features  global 
optimization  with  real-space  representation  and  flexible  physically  motivated  variation  operators. 
For  each  candidate  structure  generated  by  USPEX,  the  density  functional  theory  (DFT) 
calculations  were  performed  at  OK  by  application  of  the  Vienna  Ab-initio  Simulation  Package 
(VASP)  [28,  29].  The  compounds  in  the  Ta-C  system  explored  in  this  work  contained  up  to  30 
atoms  in  the  unit  cell.  The  initial  generation  consisted  of  50  structures  produced  randomly.  All  the 
subsequent  generations  contained  40  structures,  of  these,  40%  of  the  new  structures  were  produced 
by  heredity,  20%  produced  by  soft  mutation,  another  20%  produced  by  transmutation,  and  the  final 
20%  produced  randomly.  Thermodynamically  stable  structures  were  detennined  by  a  convex  hull 
construction;  a  phase  was  stable  if  its  decomposition  enthalpy  into  any  other  structures  was 
positive  [30].  The  calculation  was  tenninated  when  the  stable  structures  did  not  change  for  10 
generations.  The  VASP  simulations  utilized  the  projector  augmented  wave  method  [31,32]  in 
conjunction  with  the  generalized  gradient  approximation  (GGA)  using  the  parameterization  by 
Perdew,  Burke  and  Ernzerhof  (PBE)  [33],  Plane  waves  with  a  cutoff  energy  of  600  eV  and  k-point 
mesh  resolution  in  reciprocal  space  of  27ix0.03  A”1  yielded  converged  results.  Denser  k-point 
meshes  in  reciprocal  space,  with  a  resolution  of  27tx0.02A~1  were  employed  for  the  elastic 
constants  and  density  of  states  (DOS)  calculations.  The  phonon  spectra  were  obtained  by  Phonopy 
code  [34]  using  the  force  constants  calculated  by  VASP. 
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Results  and  discussions 

Thermal  stability  of  various  compounds 

The  enthalpies  of  formation  for  each  of  the  compounds  identified  are  shown  in  figure  2(a), 
while  the  stable  compounds  are  highlighted  in  black  and  tabulated  in  Table  I  and  II.  The 
thermodynamically  stable  phases  were  found  to  be  TaC  (Fm3m),  TaeCs  (C2/m)  and  TaiC  (P3ml). 
Since  TaC  and  Ta2C  are  commonly  observed  in  the  literature,  these  results  were  not  unexpected. 
In  contrast,  the  stability  of  TaeCs  could  be  considered,  by  some,  a  surprise.  No  imaginary  phonon 
frequencies  were  found  throughout  TaeCs’ s  Brillouin  zone’s  phonon  dispersion  curves,  figure  2, 
suggesting  dynamical  stability.  For  comparison,  the  formation  enthalpies  of  the  two  proposed 
TaeCs  structures  (C2  and  P3i)  by  Gusev  [14]  were  also  calculated  and  tabulated  in  Table  I  and  II. 
Their  phonon  spectrums  also  were  found  not  to  have  any  imaginary  frequency  and  are  also  plotted 
figure  3.  Interestingly,  these  results  demonstrate  that  these  TaeCs  phases  cannot  be  distinguished 
energetically  by  DFT  calculations,  which  means  they  are  degenerate.  The  consequences  of  this 
will  be  further  developed  in  the  proceeding  discussion.  The  prior  C2  and  P3i  structures  were 
assumed  to  have  vacancies  in  perfect  sites;  however,  after  structural  relaxation,  the  structures 
changed  slightly.  The  C2  structure  became  C2/c  (15)  and  the  P3ichange  to  P3 1 12  (151). 
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FIGURE  2.  (a)  The  convex  hull  construction  of  the  Ta-C  system  at  ambient  pressure.  Circles 
denote  different  structures;  those  filled  circles  located  on  the  convex  hull  are  thermodynamically 
stable  and  are  labeled  in  the  image;  the  open  circles  represent  structures  with  higher  energy.  The 
dashed  arrow  in  the  enlarged  image  points  to  the  fault-formed  TaeCs  phase  which  is  at  a  slightly 
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higher  energy  than  the  vacancy  formed  T aeCs  filled  circle  (b)  The  convex  hull  construction  of  the 
Ta-C  system  without  TaeCs.  The  Ta4C3  phase,  commonly  observed  in  experiments,  being 
identified  as  the  next  lowest  energy  structure  for  substoichiometric  TaCx  to  decompose  into.  The 
dashed  arrow  in  the  enlarged  image  points  to  a  vacancy  ordered  Ta4C3  phase  that  is  at  a  slightly 
higher  energy  than  the  fault-fonned  Ta4C3  phase  that  is  now  in  contact  with  the  convex  hull 
boundary. 


Wave  vector 


FIGURE  3.  Phonon  dispersion  curves  of  TaeCs:  (a)  C2/m  (b)  C2/c  (c)  P3 1 12. 
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Table  I.  The  crystallographic  data  of  stable  and  meta-stable  tantalum  carbides. 


Compound 

Space 

group 

This  work 

Lattice  constants  (A) 

Experimental  results 

Atom  positions 
(Wyckoff  position) 

TaC 

Fm3m 

a=4.47 

a=4.46[41] 

Ta(4a)  (0.0,0.0,0.0) 

C(4b)  (0.5, 0.5, 0.5) 

TaeCs 

C2/m 

a=  5.46 
b=  9.47 
c=  5.45 
|3=109.67° 

Ta(8j)  (0.24358,0.17417,0.24757) 

Ta(4i)  (0.76074,0.00000,0.23907) 

C(4g)  (0.00000,0.66677,0.00000) 

C(4h)  (0.00000,0.16436,0.50000) 

C(2a)  (0.00000,0.00000,0.00000) 

C2/c 

a=  5.47 
b=  9.47 
c=  10.89 
P=109.41° 

Ta(8f)  (0.75926,  0.00744,  0.8758) 

Ta(8f)(0 .735 14,  0.33269,  0.8761) 

Ta(8f)  (0.26227,  0.15857 , 0.8814) 

C(4e)  (0.0,  0.33354 , 0.25) 

C(4e)  (0.0,  0.998  ,  0.25) 

C(4b)  (0.0,  0.5 , 0.0) 

C(8f)  (0.5,  0.333 , 0.0) 

P3 1 12 

a=  5.46 
b=  5.46 
c=  15.42 
y=120.0o 

Ta(6c)  (0.1 1207,0.20788,  0.41580  ) 

Ta(6c)  (0.45943,  -0.0944,  0.41253) 

Ta(6c)  (0.76274,  0.55730,0.41613) 

C(3b)  (0.77775,  0.22225,  0.83333) 

C(3b)  (0.1 1333,  0.88667,  0.83333  ) 

C(3a)  (0.1 1062,  0.88938,  0.33333) 

C(3a)  (0.44351,  0.55649,  0.33333  ) 

C(3a)  (0.77745,  0.22255,  0.33333) 

R3m 

a=3.15 
b=  3.15 
c=45.84 
7=120.0° 

Ta(6c)  (0.000,  0.000,0.58530) 

Ta(6c)  (0.000,  0.000,0.30647) 

Ta(6c)  (0.000,  0.000,0.13907) 

C(6c)  (0.000,  0.000,  -0.05305) 

C(6c)  (0.000,  0.000,  0.77664) 

C(3b)  (0.000,  0.000,  0.50000) 

TatCs 

R3m 

a=3.14 
b=  3.14 
c=  30.39 
7=120.0° 

a=3.12,c=  30.06  [19] 

Ta(6c)  (0.000,  0.000,  0.78952) 

Ta(6c)  (0.000,  0.000,  0.62592) 

C(3b)  (0.000,  0.000,  0.50000) 

C(6c)  (0.000,  0.000,  0.08036) 

C2/c 

a=  6.23 
b=  12.56 
c=  5.45 
P=124.86° 

Ta(8f)  (0.24822,0.19059,0.75122) 

Ta(8f)  (0.23955,  0.06458,  0.23260) 

C(4e)  (0.000,  0.18464,  0.2500) 

C(4e)  (0.000,  -0.06186,  0.2500) 

C(4e)  (0.000,  0.68900,  0.2500) 

Ta3C2 

R3m 

a=3.14 
b=  3.14 
c=  22.45 
7=120.0° 

Ta(3a)  (0.000,  0.000,  0.00000) 

Ta(6c)  (0.000,  0.000,  0.22069) 

C(6c)  (0.000,  0.000,  0.60798) 

Ta2C 

P3ml 

a=3.1 1, 
b=3.11 
c=  4.95 
7=120.0° 

a=3.10,  c=  4.94[42] 

Ta(12j)  (0.333,0.667,0.753) 

C(6h)  (0.0, 0.0, 0.5) 
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Table  II.  The  formation  enthalpies  of  the  stable  and  meta-stable  tantalum  carbides. 


Compound 

Space  group 

Enthalpy  of  fonnation  (eV/atom) 

TaC 

Fm3m 

-0.59 

TaeCs 

C2/m 

-0.63 

C2/c 

-0.63 

P3  i  12 

-0.63 

R3m 

-0.59 

TatCs 

R3m 

-0.62 

C2/c 

-0.61 

Ta3C2 

R3m 

-0.61 

Ta2C 

P3ml 

-0.60 

As  mentioned  by  Gusev  [14],  these  superstructures  have  identical  short-range  order  and  are 
similar  up  to  the  third  coordination  sphere  with  respect  to  the  metal  atom.  Consequently,  they  have 
nearly  equal  fonnation  energies  as  proven  by  this  work.  Unlike  the  niobium  and  vanadium 
carbides,  where  the  MeeCs  phase  is  observed  more  frequently  in  experiments  [18],  the  lack  of 
consistent  identification  of  this  phase  in  the  tantalum  carbides  could  be  attributed  to  several 
(potentially  compounding)  factors.  First,  as  noted  by  Gusev  [19]  and  alluded  to  above,  this  phase 
appears  to  require  very  slow  cooling  rates  in  order  for  vacancy  ordering  to  occur.  These  cooling 
rates  are  impractical  in  most  high  temperature  fabrication  processing  resulting  in  the  higher 
temperature  rocksalt  TaCi-x  to  be  quenched-in  at  room  temperature.  Second,  the  vacancy  migration 
energies  (formation  +  diffusion)  within  the  tantalum  carbides  (-4.3  eV)  are  higher  than  either  the 
niobium  carbides  (~3.6  eV)  and  vanadium  carbides  (-3.0  eV)  [35],  where  these  phases  are  more 
readily  observed.  These  computed  vacancy  formation  energies  were  found  from  the  energy  needed 
to  extract  a  carbon  atom  from  the  B 1  metal  carbide  structure;  the  migration  energy  was  computed 
from  the  minimum  energy  path  determined  using  the  nudged-elastic-band  (NEB)  method.  This 
higher  energy  barrier  would  potentially  further  hinder  phase  transformations,  which  occurs 
through  long-range  vacancy  ordering.  Supporting  this  idea  is  the  experimental  observations  in 
niobium  and  vanadium  carbides  that  MeeCs  fonns  as  domains  within  the  monocarbide  rather  than 
consuming  the  full  grain  [36].  Third,  the  indistinguishable  enthalpies  of  fonnation  of  the  TaeCs 
structures,  Table  II,  suggest  that  multiple  types  of  long  range  ordering  may  occur,  which  would 
make  it  difficult  to  identify  a  specific  TaeCs  structure.  Finally,  the  larger  scattering  ratio  of 
tantalum-to-carbon,  in  comparison  to  other  group  VB  transition  metal  carbides,  and  similar 
electron  diffraction  patterns  further  inhibits  clear  identification  and  classification  of  the  MeeCs 
phase  in  the  tantalum  carbon  system. 

The  second  lowest  energy  TaeCs  structure  found  by  USPEX  is  the  R3m  phase,  also  tabulated 
in  Table  I  and  II.  It’s  location  on  the  convex  hull  plot  of  figure  2(a)  is  indicated  by  the  dashed 
arrow.  This  is  a  ‘faulted  structure’  similar  to  the  Ta4C3  phase  that  will  be  further  discussed  below. 
The  formation  enthalpy  of  this  phase  is  AE  -0.04  eV/atom  higher  than  those  of  the  vacancy  ordered 
fonns.  Considering  that  this  calculation  was  done  at  OK,  this  small  energy  difference  can  be  easily 
overcome  with  increasing  temperature  by  the  vibration  contributions  to  the  free  energy  of 
fonnation.  Though  vacancy  ordered  phases  have  slightly  lower  fonnation  enthalpies,  the  kinetic 
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process  of  ordering  also  makes  them  possibly  more  difficult  to  occur.  In  contrast,  the  faulted 
structure  is  potentially  more  likely  because  the  fault  formation  process  requires  vacancy 
condensation  to  one  particular  plane,  a  localized  effect,  which  would  be  followed  by  shear 
associated  with  Shockley  dislocations  to  create  the  fault  [13].  The  true  formation  of  this  faulted 
TaeCs  phase  would  then  require  the  appropriate  number  and  spacing  of  these  types  of  faults.  In 
addition,  the  competition  between  the  ‘faulted’  structure  and  the  vacancy  ordered  forms  is  further 
complicated  by  kinetic  processes. 

Besides  the  aforementioned  stable  compounds,  Ta4C3  (R3m)  and  Ta3C2(R3m)  were  found  to 
be  very  close  to  the  convex  hull  curve,  figure  2(a),  and  whose  enthalpies  of  formation  are  slightly 
higher  than  the  thermodynamic  stable  phases,  Table  II.  Recognizing  that  the  TaeCs  phase  is  not 
commonly  observed  experimentally,  it  has  been  removed  and  a  ‘metastable’  convex  hull  boundary 
is  re-plotted  in  figure  2(b).  In  this  case,  the  Ta4C3,  along  with  TaC  and  Ta2C,  are  on  the  convex 
hull  with  Ta3C2  at  a  slightly  higher  energy.  The  slightly  higher  energy  of  Ta3C2may  explain  why 
this  phase,  as  compared  to  Ta4C3,  has  not  been  reported  experimentally. 

These  metal-rich  carbide  fault-formed  structures  (R3m)  can  be  created  from  the  B1  TaC 
structure  by  the  removal  of  specific  {111}  carbon  layers.  This  carbon  removal  would  then  place 
metal-metal  Ta  {111}  planes  into  contact.  A  subsequent  shear  on  this  metal-metal  contact  plane 
created  by  Shockley  partial  dislocations  would  then  create  the  requisite  stacking  fault.  The 
conventional  trigonal  unit  cell  for  each  of  these  phases  must  contain  3  faults  to  match  the  stacking 
sequence  of  the  parent  TaC  phase.  This  results  in  there  being  3  times  as  many  metal  planes  in  the 
stacking  sequence  as  are  required  for  stoichiometry  for  this  phase.  Specifically,  there  are  12  metal 
layers  in  Ta4C3,  9  in  Ta3C2  and  18  in  TaeCs  as  shown  in  Table  III.  The  resulting  stacking  sequence 
of  Ta3C2  isAyBaC|BaC|3A|C[3AyB  where  A,B  and  C  denote  the  stacking  of  the  Ta 
atoms;  a,  P  and  y  denote  the  stacking  of  the  carbon  atoms;  and  |  denotes  a  removal  of  a  carbon 
plane  and  creation  of  the  stacking  fault.  This  fault-based  formation  mechanism  proposed  by 
Rowcliff  [13]  provides  the  removal  of  the  carbon  for  the  local  composition  change  with  the  shear 
necessary  for  the  crystallographic  change  that  creates  the  new  phase.  It  is  worth  pointing  out  that 
Ta2C  follows  this  carbon  removal  strategy  but  only  has  two  metal  atoms  in  the  trigonal  unit  cell 
and  belongs  to  a  different  space  group,  P3ml;  this  is  because  of  its  unique  stacking.  The  model 
given  above  collapses  into  the  higher  symmetry  P3ml  class  though  they  share  the  same  point 
group  3. 

The  ordered  removal  of  a  carbon  plane,  or  alternatively  the  ordering  of  vacancies  on  a  plane, 
is  similar  to  experimental  observations  of  Frank  loops  in  face  centered  cubic  metals  [37]  or  the 
precipitation  mechanisms  of  y’  in  a  in  Al-Ag  alloys  [38].  Wang  et  al.  [39]  recently  reported  that 
the  generalized  stacking  fault  energy  in  Ta2C  is  an  order  of  magnitude  lower  for  a  Ta-Ta  verses  a 
Ta-C  close-packed  plane  arrangement.  This  result  would  support  that  the  energetically  favorable 
shear  would  follow  the  removal  of  the  carbon  plane.  For  the  faulted  TaeCs  phase  to  form,  every 
sixth  close  packed  carbon  plane  would  be  removed  and  sheared  in  sequential  sequence  with  every 
five  Ta-C  stacks  to  create  the  unit  cell;  similarly,  every  fourth  close  packed  carbon  plane  would  be 
removed  and  sheared  in  sequential  sequence  with  every  three  Ta-C  stacks  to  create  the  unit  cell  of 
Ta4C3;  and  so  forth.  This  removal  and  shear  patterning  provides  qualitative  insights  into  why  the 
TaeCs  phase  (or  higher  ordered  fault-forming  phases)  may  not  readily  form  as  it  requires  longer- 
range  vacancy  ordering  over  multiple  sequences  of  Ta/C  layered  planes. 
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Considering  the  enthalpies  of  fonnation  of  the  structures  found  from  the  USPEX  codes,  it  is 
clear  there  is  a  competition  between  the  vacancy  ordered  forms  of  the  B 1  structure  and  those  phases 
that  are  essentially  formed  from  stacking  faults.  At  the  composition  of  TaeCsThe  faulted  structures 
is  0.04  eV  higher  than  the  vacancy  ordered  structure.  However,  at  the  Ta4C3  composition,  the 
faulted  zeta-phase  fonns  as  described  above.  In  examining  the  DFT  results,  it  is  found  that  the 
second  lowest  energy  Ta4C3  structure  is  the  C2/c  vacancy  ordered  phase,  Table  II,  and  is  indicated 
by  its  location  on  the  ‘metastable’  convex  hull  plot  by  the  dashed  arrow  in  figure  2(b).  This  clearly 
indicates  an  energy  competition  present  in  these  substoichiometric  compounds  between  fault¬ 
forming  and  vacancy  order  rock-salt  based  phases,  with  the  fault-forming  structures  becoming 
more  thermodynamically  favorable  with  further  loss  of  carbon. 

Table  III.  The  stacking  sequences  of  all  the  faulted  structures  in  tantalum  carbides.  The  Roman 
letters  represent  the  metal  close-pack  stacking  sequence  and  the  Greek  letters  represent  the 

equivalent  for  the  carbon  layers. 


TaC 

TaeCs 

T34C3 

Ta3C2 

Ta2C 

A 

A 

A 

A 

A 

y 

y 

y 

y 

y 

B 

B 

B 

B 

B 

a 

Stacking  fault 

Stacking  fault 

Stacking  fault 

C 

A 

A 

A 

P 

y 

y 
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A 

B 

B 

B 
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C 

P 

P 
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The  new  ‘metastable’  convex  hull  construction  in  figure  2(b)  also  sheds  light  on  prior 
experimental  data  where  TaC,  Ta4C3,  and  TaiC  were  all  found  to  co-exist  in  a  63Ta:37C  (at.  %) 
compound,  figure  fib)  [1].  The  energies  of  each  of  these  systems  are  very  near  each  other  and  all 
are  at  the  minimum  when  the  kinetically  constrained  TaeCs  vacancy  ordered  phase  does  not  fonn. 
Since  each  of  these  phases  are  related  to  each  other  through  the  loss  of  carbon  and  formation  of 
stacking  faults,  it  becomes  very  difficult  to  distinguish  between  a  true  Ta4C3  phase  and  a  rocksalt 
phase  with  a  set  of  random  stacking  faults.  Consequently,  the  decomposition  of  the  TaCi-xis  not 
strongly  driven  to  one  particular  phase.  As  alluded  to  above,  the  collective  system  is  ‘frustrated’ 
in  which  phase  to  form.  This  idea  is  supported  partially  by  experimental  observations  of  the  phase 
morphology  near  the  Ta4C3  composition  where  the  material  is  observed  to  be  composed  of  very 
fine  planar  morphologies.  Though  these  planar  laths  may  appear  thick  in  two-dimensional 
projection  micrographs,  figure  1(a),  Morris  et  al.  provided  transmission  electron  tomography 
revealing  their  fine  thickness  and  plate  like  structure  in  the  plane  of  the  carbide  matrix  [  1] .  Detailed 
microscopy  work  in  TaiC  [40]  has  also  demonstrated  that  slightly  carbon-rich  Ta2C  will  also  create 
a  series  of  faults  whose  stacking  sequence  agrees  with  the  phase  transformation  from  TaiC  through 
Ta4C3  to  TaC.  The  phase  morphology  and  phase  transformation  pathways  observed  experimentally 
are  very  likely  promoted  by  the  close  energetics  of  the  TaiC,  Ta3C2  and  Ta4C3  phases  revealed 
through  the  computed  work  here. 

The  evolutionary  algorithm  has  been  able  to  correctly  identify  the  experimentally  determined 
phases  as  well  as  the  metastable  phases.  Considering  the  small  energy  differences  between  several 
of  these  phases,  it  is  somewhat  remarkable  to  consider  how  well  the  thermodynamic  stability  of 
these  high-temperature  materials  is  described  by  the  used  T=0K  method.  Regardless,  this  has 
provided  new  understanding  between  the  thennodynamic  competition  of  the  vacancy  ordered  and 
stacking  fault  variations  of  these  carbide  phases  that  results  in  the  complex  microstructures 
reported. 

Elastic  properties 

As  noted  in  the  introduction  section,  the  tantalum  carbides  exhibit  a  range  of  mechanical 
responses.  Through  tailoring  the  phase  content  and  microstructure,  these  attributes  can  be  tuned 
for  various  loading  demands.  Since  the  production  of  single  crystals  can  be  arduous,  if  not 
impossible,  for  some  of  the  compounds,  computational  calculations  can  provide  an  alternative  way 
to  probe  the  elastic  properties  of  these  phases.  The  elastic  constants  for  each  of  the  stable  and 
metastable  compounds  have  been  computed,  and  compared  where  possible  to  either  experimental 
and/or  other  computational  investigations  [15,43].  The  tabulation  of  the  elastic  constants  is  found 
in  Table  IV.  For  cubic  TaC,  there  are  3  nonzero  independent  elastic  constants;  for  monoclinic 
C2/m  and  C2/c  TaeCs,  there  are  13;  for  P3 1 12  and  P3ml  with  the  trigonal  symmetry,  there  are  7 
and  6  respectively  [44]. 
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TABLE  IV.  Calculated  elastic  constants  C,h  the  bulk  modulus,  shear  modulus  and  hardness 

(GPa)  ofTaC,Ta6C5,  Ta4C3,  Ta3C2,and  Ta2C. 


Compound 

TaC 

Ta6C5 

Ta4C3 

Ta3C2 

Ta2C 

Space  group 

Fm3m 

C2/m 

C2/c 

P3i  12 

R3m 

R3m 

C2/c 

R3m 

P3ml 

Cn(GPa) 

737 

494 

497 

499 

585 

571 

504 

521 

479 

C22(GPa) 

495 

500 

493 

C33(GPa) 

479 

487 

486 

549 

537 

540 

504 

500 

C44(GPa) 

175 

173 

169 

206 

211 

179 

179 

164 

133 

C55(GPa) 

202 

209 

138 

C66(GPa) 

213 

206 

140 

Ci2(GPa) 

141 

157 

146 

150 

164 

156 

202 

138 

164 

CnCGPa) 

177 

178 

179 

188 

169 

130 

192 

149 

Ci4(GPa) 

-27 

-49 

47 

-44 

-45 

Ci5(GPa) 

Ci6(GPa) 

-38 

-38 

27 

-0.8 

C23(GPa) 

192 

178 

120 

C26(GPa) 

34 

36 

-0.7 

C36(GPa) 

5 

-2 

11.4 

C45(GPa) 

41 

35 

14.2 

B(GPa) 

339.67 

279.96 

276.39 

277.75 

311 

296 

271 

288 

264.66 

G(GPa) 

216.90 

175.87 

178.40 

179.71 

200.4 

187 

161 

167 

148.08 

k 

0.64 

0.63 

0.65 

0.65 

0.64 

0.63 

0.59 

0.58 

0.56 

Hv(GPa) 

24.53 

20.89 

21.87 

22.05 

23.6 

21.92 

18.3 

18.19 

15.87 

In  order  to  easily  compare  the  trends  in  hardness  among  the  carbides,  we  estimated  the  elastic 
polycrystalline  elastic  constants  using  the  well-known  Voigt  ( Gv )  and  Reuss  (Gr)  elastic  constants 
method  [45,46,47]  with  the  average  shear  modulus,  G,  defined  as 

G  =  ~(GV  +  Gr)[45]. 

The  “theoretical”  hardness,  Hv,  for  all  the  carbides  was  then  estimated  using 

Hv  =  2(k2G)0,585-3  [48] 

where  k  is  Pugh's  modulus  ratio  defined  as  G/B  with  B  being  the  bulk  modulus.  The  hardness  of 
the  Ta-C  phase  as  a  function  of  carbon  content  is  plotted  in  figure  4.  This  is  referred  as  “theoretical” 
hardness  because  it  is  based  solely  on  the  computed  elastic  constants  and  lacks  any  contribution 
of  microstructure  to  hardness.  It  can  be  seen  that  all  of  the  carbides  elastically  soften  and  the 
hardness  monotonically  decreases  with  decreasing  carbon  content. 

Interestingly,  one  thing  that  that  this  derived  hardness  does  not  predict  is  the  anomalous 
increase  in  hardness  reported  for  substoichiometric  compositions  in  tantalum,  niobium,  and 
vanadium  carbides  [49,50].  This  increase  in  hardness  appears  near  the  MeeCs  composition;  hence 
the  elastic  constants  are  insufficient  to  fully  elucidate  the  experimental  hardness  behavior,  a  point 
that  we  demonstrate  concretely  in  a  later  work  [35].  These  results  demonstrate  that  the  elastic 
constants  of  ordered  tantalum  carbides  decrease  continually  with  carbon  loss,  which  corresponds 
to  a  decrease  in  covalent  bonding  which  is  demonstrated  below. 
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FIGURE  4.  The  theoretical  hardness  for  the  tantalum  carbides  calculated  in  this  paper. 
Electronic  structure 


In  a  prior  computational  study  of  TaiC  [39],  the  bonding  structure  was  found  to  be 
paramount  in  dictating  the  slip  behavior.  In  this  particular  system,  a  local  metallic  bond  was  present 
and  allowed,  in  specific  <hkl>  { hkl }  orientations,  an  order  of  magnitude  decreases  in  the 
generalized  stacking  fault  energy  which  contributed  to  the  enhanced  plasticity.  Thus,  due  to  the 
correlation  between  bonding  and  properties,  it  is  important  to  characterize  the  bonding  in  each  of 
the  stable  tantalum  carbide  structures.  Using  the  electronic  DOS,  the  bonding  characteristics  in 
each  thermodynamic  stable  tantalum-carbide  was  investigated  and  plotted  in  figure  5. 


Energy  (eV) 


FIGURE  5.  The  total  density  of  states  for  TaC,  TaeCsand  Ta2Ci.  The  DOS  is  normalized  by  the 
number  of  unit  formula  units  in  the  simulation  and  the  Fermi  level  has  been  shifted  to  0  eV. 
(Color  available  on  line) 
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The  peak  in  the  DOS  for  TaC  is  at  -  -5  eV  and  is  associated  with  the  pd  bonding  between 
the  tantalum  and  carbon  atoms,  which  is  covalent  in  nature  [43],  It  appears  that  this  peak  decreases 
with  carbon  content,  which  implies  it  has  less  contribution  to  the  bonding  characteristics  in  the 
metal-rich  carbides.  Similarly,  the  density  of  the  d  state  increases  near  the  Fenni  level,  which  can 
be  interpreted  as  an  increase  in  the  dd  bonding  or  metallic  nature  bonds  between  the  tantalum 
atoms.  This  is  the  highest  for  the  Ta2C  system.  The  change  between  the p  electrons  (covalent)  and 
the  d  electrons  (metallic)  with  carbon  content  results  in  a  collective  decrease  in  bond  strength  for 
the  carbides.  This  corresponds  to  the  decrease  in  elastic  constants  tabulated  for  the  metal-rich 
carbides  in  Table  III.  In  general,  the  melting  temperature  of  a  material  scales  with  the  bond 
strength.  This  trend  is  also  evident  in  this  plot.  TaC  melts  at  -4200 K  where  as  Ta2C  melts  at 
3600K,  with  TaC  having  more  of  a  covalent  bond  character.  The  TaeCs  phase,  according  to  Gusev 
[19],  is  lost  at  -143  OK,  which  would  deviate  from  this  trend,  expect  that  this  phase  is  a  vacancy 
ordered  phase  and  does  not  melt  but  rather  becomes  disordered.  The  computational  findings  appear 
to  trend  well  with  the  physical  properties  of  these  materials. 

Conclusions 

By  using  a  variable-composition  ab  initio  evolutionary  algorithm,  we  explored  the  stable  and 
metastable  compounds  in  the  Ta-C  system  at  ambient  pressure  and  OK.  The  results  revealed  that 
TaC,  TaeCs  and  Ta2C  are  the  thermodynamically  lowest  energy  configurations  with  Ta4C3  and 
Ta3C2  having  slightly  higher  energies.  The  vacancy  ordered  fonn  of  TaeCs  had  three  degenerate 
energetic  structures  (C2/m,  C2/c,  and  P3 1 12).  In  addition,  the  stacking-fault  based  TaeCs  structure 
has  a  slightly  higher  energy,  AE-0.04  eV/atom  than  the  vacancy  ordered  forms  with  an  equivalent 
space  group,  R3m,  to  the  Ta4C3  phase.  Such  a  small  energy  difference  could  be  overcome  with 
increasing  thennal  contributions,  which  are  not  considered  in  this  computational  algorithm 
approach.  Regardless,  the  computational  methodology  provided  thermodynamic  stability  energies 
that  are  agreement  with  the  phases  noted  in  these  high-temperature  materials.  In  contrast  to  TaeCs, 
the  Ta4C3  phase  is  commonly  observed  to  form  in  experiments.  By  removing  the  kinetically 
constrained  TaeCs  phase  from  the  convex  hull  plot,  it  was  re -plotted  revealing  that  TaC,  Ta4C3, 
and  Ta2C  were  the  lowest  ‘stable’  phases,  which  are  all  commonly  observed.  The  evolutional 
algorithm  was  successful  at  predicting  all  of  these  experimentally  detennined  phases.  Interestingly, 
these  simulations  demonstrated  that  in  addition  to  the  fault-formed  R3m  Ta4C3  phase,  there  exists 
a  slightly  higher  energy,  AE-0.01  eV/atom,  vacancy  ordered  C2/c  phase.  From  these 
computational  simulations,  a  thermodynamic  competition  was  noted  to  occur  between  the  vacancy 
ordering  and  fault-fonned  phases,  with  the  latter  structures  being  more  favorable  with  decreasing 
carbon  content  in  TaCi-x.  While  the  stable  phase  predicted  here  are  in  agreement  with  the  most 
comprehensive  phase  diagrams  of  Ta-C  by  Gusev  [19],  these  simulations  provide  additional 
insight  into  the  competition  between  the  vacancy  ordered  B 1  phases  and  fault  forming  phases  that 
are  not  readily  apparent  from  experiments. 

Collectively,  the  close  enthalpy  of  fonnation  energies  computed  from  this  study  for  each  of 
the  stacking-fault  based  carbides  may  ‘frustrate’  the  carbide  in  the  co-precipitation  of  multiple 
phases  in  substoichiometric  tantalum  carbides  which  leads  to  the  complex  microstructures. 

The  DFT  calculations  also  provided  the  elastic  constants  for  each  of  the  stable  and  metastable 
phases.  As  the  carbon  content  increased,  the  elastic  constants  values  increased.  This  has  been 
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associated  with  the  increase  in  the  pd  covalent  and  decrease  in  the  metallic  dd  bonding  of  the 
carbide  from  the  density  of  states. 
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Dissemination  of  Research  Findings: 

MMM  2014  -  7th  International  Conference  on  Multiscale  Materials  Modeling  conference,  Oct. 

6- 1 0th,  2014,  Berkeley,  CA;  mmm2014berkeley.iop.org 
Oral  presentation : 

>  “Metal-rich  Ceramic  Phase  Stability  in  Group  IV  and  V  Carbides  and  Nitrides” 
Xiao-Xiang  Yu,  Christopher  R.  Weinberger,  Gregory  B.  Thompson 

>  “Linking  Computational  Stacking  Fault  Energies  with  Plasticity  Mechanisms  in  High 
Temperature  Carbide  Ceramics” 

Gregory  B.  Thompson,  Nicholas  De  Leon,  Xiao-Xiang  Yu,  Billie  Wang,  and  Christopher 
R.  Weinberger 

Poster  presentation : 

>  “Computational  Investigations  into  the  Hardness  of  Metal  Carbides” 

Xiao-Xiang  Yu,  Christopher  R.  Weinberger,  Gregory  B.  Thompson 

Materials  Science  and  Technology  (MS&T)  Oct.  12-16,  2014,  Pittsburg,  PA  matscitech.org 

Oral  presentation 

r  “Ab  Initio  Investigations  of  the  Phase  Stability  and  Hardness  in  Tantalum  Carbides” 
Xiao-Xiang  Yu,  Christopher  R.  Weinberger,  Gregory  B.  Thompson 

United  States  Advanced  Ceramics  Association  (USACA)  January  26th  -  30th,  2015,  Coco 
Beach,  FL 

Oral  presentation 

>  “Phase  Stability  in  Tantalum  Carbides  and  its  Relationship  to  Microstructure  and 
Mechanical  Behavior”  Xiao-Xiang  Yu,  Christopher  R.  Weinberger,  and  Gregory  B. 
Thompson 

The  American  Ceramic  Society  International  Conference  and  Expo  on  Advanced  Ceramics 
and  Composites  January  25th  -  30th,  2015,  Daytona,  FL 

Oral  presentation 

>  “Understanding  Phase  Stability  in  the  Transition  Metal  Group  IVB  and  VB  Carbides” 
Xiao-Xiang  Yu,  Christopher  R.  Weinberger,  and  Gregory  B.  Thompson 
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The  Minerals,  Metals,  and  Materials  Society  (TMS)  144th  Annual  Meeting  &  Exhbinition, 
March  15-1 9th  2015,  Orlando,  FL 

Oral  presentation 


>  “First-Principles  Investigation  on  the  Phase  Stability  in  Tantalum  Carbides”  Xiao-Xiang 
Yu,  Christopher  R.  Weinberger,  and  Gregory  B.  Thompson 

Poster  presentation : 

>  “A  Non-Contact  Means  of  Thermo-Mechanical  Testing  of  Materials  beyond  2000  °C” 
Chase  Smith,  Nicholas  De  Feon,  and  Gregory  Thompson 


Honors: 


PI  Thompson  was  the  recipient  of  the  2014  Blackmon-Moody  Outstanding  Professor  of  the  Year 
award  for  UA  on  October  24th,  2014.  Below  is  a  picture  from  the  award  ceremony.  More  info 
can  be  found  below  or  with  the  proceeding  article. 


http://uanews.ua.edu/2014/10/engineering-professor-to-receive-uas-blackmon-moody-award/ 


Mr.  Blackmon  (left);  President  of  UA,  Dr.  Judy  Bonner  (center);  Professor  Greg  Thompson 
(right)  receiving  the  award 
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Engineering  Professor  to  Receive  UA’s  Blackmon-Moody  Award 

TUSCALOOSA,  Ala.  —  Dr.  Gregory  B.  Thompson,  professor  of  metallurgical  and  materials 
engineering  at  The  University  of  Alabama,  will  receive  the  2014  Blackmon-Moody  Outstanding 
Professor  Award  at  an  Oct.  24  ceremony. 

Since  coming  to  the  University  in  2003,  Thompson  has  received  external  funding  from 
organizations  such  as  the  National  Science  Foundation,  NASA,  the  Air  Force  Office  of  Scientific 
Research,  the  Office  of  Naval  Research  and  the  Army  Research  Office.  These  awards  have 
totaled  several  millions  of  dollars  in  the  support  of  basic  and  applied  research  in  materials 
development. 

Thompson  serves  as  director  of  the  Materials  Science  doctoral  program  on  campus,  associate 
director  of  the  University’s  Central  Analytical  Facility  and  director  of  the  Materials  Metrology 
Research  Consortium. 

Thompson’s  research  emphasis  is  in  the  thermodynamics  and  mechanisms  of  phase 
transformations  and  their  effect  on  microstructure.  In  recent  years,  his  research  group  has  studied 
different  classes  of  transition  metal  ceramics  and  the  levels  of  ductility  they  achieve  at  elevated 
temperatures. 

His  research  uses  a  variety  of  state-of-the-art  analytical  microscopes  that  characterize  materials 
from  the  atomic  to  micron-length  scales.  He  has  written  more  than  100  peer-reviewed  articles 
and  three  book  chapters  in  his  research  areas  of  analytical  microscopy  and  phase  transformations. 

Thompson  has  been  involved  in  multiple  outreach  and  service  activities.  In  2013  and  2014,  he 
was  a  co-organizer  for  an  ASM  Materials  Camp  on  campus,  a  program  aimed  at  teaching 
middle-  and  high-school  teachers  how  to  incorporate  materials  engineering  into  secondary 
education  classrooms. 

In  2004,  he  co-founded  the  Nanoscience  and  Engineering  High  School  Internship  Program,  a  10- 
week  summer  program  in  which  high-school  students  work  with  UA  faculty  on  topics  like 
magnetic  nanoparticles,  sensor  technology  and  catalysis. 

In  2012,  Thompson  served  as  organizing  chair  for  the  53rd  International  Field  Emission 
Symposium  that  brought  more  than  170  scientists  from  15  countries  to  the  campus  for  a  week- 
long  conference  on  atom  probe  and  high  field  nanosciences. 

In  2013,  Thompson  was  named  a  Crimson  Tide  Hometown  Hero  for  his  involvement  in  relief 
efforts  following  the  tornado  that  struck  Tuscaloosa  April  27,  2011.  Additionally,  he  was  been 
recognized  multiple  times  by  the  Minerals,  Metals  and  Materials  Society,  having  been  honored 
as  a  Young  Leader  in  2005  and  a  Young  Leader  International  Scholar  in  2008. 

Earning  his  bachelor’s  degree  in  physics  from  Brigham  Young  University  in  Provo,  Utah,  in 
1996,  he  obtained  both  his  master’s  degree  and  his  doctorate  in  materials  science  and 
engineering  from  the  Ohio  State  University  in  1998  and  2003,  respectively. 


50 


The  Frederick  Moody  Blackmon  and  Sarah  McCorkle  Moody  Outstanding  Professor  Award  is 
presented  annually  to  a  UA  faculty  member  judged  to  have  made  extraordinary  research 
contributions  that  reflect  credit  on  the  individual,  his  or  her  field  of  study  and  on  the  University. 
It  was  created  by  Frederick  Moody  Blackmon  of  Montgomery  to  honor  the  memory  of  his 
grandmother,  Sarah  McCorkle  Moody  of  Tuscaloosa. 

The  University  of  Alabama,  a  student-centered  research  university,  is  experiencing  significant 
growth  in  both  enrollment  and  academic  quality.  This  growth,  which  is  positively  impacting  the 
campus  and  the  state's  economy,  is  in  keeping  with  UA’s  vision  to  be  the  university  of  choice  for 
the  best  and  brightest  students.  UA,  the  state's  flagship  university,  is  an  academic  community 
united  in  its  commitment  to  enhancing  the  quality  of  life  for  all  Alabamians. 
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Abstract 

Transition  metal  carbides  comprise  a  class  of  high  and  ultrahigh  melting  temperature 
materials  with  tremendous  thermo-mechanical  property  potential.  This  research 
elucidates  their  phase  stability  and  mechanical  deformation  mechanisms.  New 
accomplishments  in  the  past  year  includes  the  following: 

(1 )  Ab  initio  investigations  into  the  phase  stability  in  Ta-C.  The  results  revealed  the  most 
thermodynamically  favorable  phases  (based  on  enthalpy  of  formation  at  OK)  are  TaC, 
Ta6C5  and  Ta2C.  Interestingly,  three  degenerate  vacancy  ordered  Ta6C5  structures 
were  noted,  one  of  which  is  the  first  ever  reported.  As  the  carbon  is  reduced  from  the 
stoichiometric  monocarbide,  the  system  appears  to  undergo  a  competition  between 
precipitating  vacancy-ordered  and  faulted-based  phases,  with  the  latter  structure 
becoming  more  favorable  with  metal  enrichment.  This  work  provided  new  insights  into 
why  specific  phases  form  and  do  not  form  within  this  carbide  [Yu  et  al.  Acta  Materialia  80 
(2014)  341  ].  Ongoing  efforts  are  underway  to  determine  the  full  range  of  phase  stability  in 


the  other  group  IVB,  as  well  as  group  VB  carbides. 


(2)  Building  upon  the  computational  work,  the  program  was  able  to  address  the 
anomalous  hardness  rise  noted  in  prior  literature  in  the  group  VB  carbides  which  is 
absent  in  group  IVB  carbides.  This  behavior  is  contributed  to  the  former  carbides  having 
a  lower  vacancy  migration  energy,  with  NbC  and  VC  having  a  near  zero  and  negative 
vacancy  formation  energy  respectively.  This  favorable  vacancy  formation  and  migration 
behavior  results  in  the  formation  of  Me6C5  domains.  It  was  proposed  that  these  domains 
provide  “Hall-Petch”  strengthening.  The  higher  vacancy  formation  energy,  prevalent  in 
the  group  IVB  carbides,  likely  hindered  them  from  easily  forming  such  phase  domains 
which  prevents  the  anomalous  sub-stoichiometric  rise  in  hardness.  Hence,  the 
competition  between  elastic  constant  softening  and  microstructure-based  effects  in  the 
group  VB  carbides  is  suspected  to  cause  the  peculiar  hardness  peak  [Yu  et  al  J  of 
European  Ceramic  Society  35  (201 5)  95]. 

(3)  The  program  resolved  the  outstanding  question  of  why  group  IVB  monocarbides  slip 
on  {1 1 0}  planes  whereas  group  VB  carbides  slip  on  {1 1 1 }  at  room  temperature.  HfC 
(group  IVB)  and  TaC  (group  VB)  were  used  as  the  case  study  materials  where  upon 
Vicker’s  indentation  was  performed.  TEM  dynamical  diffraction  quantified  the  slip  systems 
with  complimentary  density  functional  theory  generalized  stacking  fault  energies  (GSFEs) 
used  to  determine  the  predicted  energy  hierarchy  for  slip.  It  was  found  that  group  VB 
carbides  can  stabilize  an  intrinsic  stacking  fault  (ISF)  on  the  {111}  plane  that  is  absent  in 
the  group  IVB  carbides.  These  results  debunk  prior  hard  sphere  based  model 
interactions  that  have  been  propagated  for  decades  to  explain  the  slip  discrepancies. 

This  ISF  allows  the  group  VB  carbides  to  circumvent  the  lower  {1 1 0}  GSEF  surface.  Also 
found  was  that  TaC,  regardless  of  GSEF  surface,  was  lower  than  HfC,  which  was 
experimentally  manifested  as  a  higher  dislocation  density.  (De  Leon  et  al.  Physical 
Review  Letters,  in  press  (2015). 

(4)  Over  the  past  year,  the  program  has  constructed  a  new  electrical  magnetic 
mechanical  testing  apparatus  that  can  provide  a  non-contact  means  of  thermo¬ 
mechanical  testing  these  carbides  up  to  their  melting  temperature.  The  improvements 
include  a  larger  field  strength  electromagnet,  which  will  enable  larger  samples  to  be 
tested.  In  addition,  a  multi-point  laser  displacement  unit  for  real-time  creep 
measurements  from  the  entire  length  of  the  bar,  not  just  the  region  of  maximum 
displacement,  has  been  added.  The  testing  apparatus  is  currently  waiting  on  a  custom 
chamber  to  house  the  testing  equipment  which  will  isolate  the  samples  from  the 
environment  preventing  deleterious  oxidation  effects  at  elevated  temperature. 

(5)  Over  the  past  year,  the  program  has  disseminated  the  findings  through  6  contributed 
talks  and  2  poster  at  various  professional  society  meetings.  The  PI  has  been  invited  to 
give  a  key  note  lecture  on  the  subject  of  carbides  at  the  forthcoming  ECI  Ultrahigh 
Temperature  Ceramics  Conference  in  April  2015. 

(6)  One  MS  student,  Chase  Smith,  graduated  under  the  program  who  will  be  continuing 
his  doctoral  studies  through  a  follow-on  effort.  One  PhD  candidate,  Nicholas  De  Leon, 
passed  his  PhD  proposal  defense  and  is  on  track  to  defend  his  dissertation  from  the  work 
collected  under  this  program  in  fall  201 5. 

(6)  The  PI  was  awarded  the  2014  Blackmon-Moody  Outstanding  Professor  of  the  Year  at 
The  University  of  Alabama.  This  award  recognizes  a  faculty  member  for  their  singular 
research  accomplishments  bringing  credit  to  the  capstone  institution. 
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